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ABSTRACT 


The  status  of  the  current  investigation  into  the  radiant  emission 
from  high  temperature  air  is  described.  A  ureatment  of  emission  by 
molecular  species  and  some  new  results  in  line  broadening  is  included. 
Results  are  presented  for  air  at  temperatures  0.5  to  5.0  eV  and 
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densities  from  1.0  to  10  times  normal. 
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I.  Introduction 

ftie  previous  report  distributed  under  this  program  (Ref.  l)  described 
our  approach  toward  a  theoretical  prediction  of  the  spectral  radiant  emission 
from  high  temperature  atomic  air  and  included  some  preliminary  results.  This 
report  describes  an  extension  of  the  previous  results  toward  lower  temperatures 
at  which  radiation  from  molecular  species  must  be  included. 

The  atomic  transitions  dominating  the  integrated  emissivity  are  gen¬ 
erally  of  reasonably  well-known  strength,  width,  and  energy  (Ref.  2  and  3)* 

Thus  we  expect  our  estimates  of  the  spectral  air  absorption  coefficient  to 
be  quite  reliable  where  the  absorption  coefficient  is  large.  Indeed,  our 
integrated  emissivities  agreed  well  with  other  reported  values  (Ref.  4,  5) 
which,  in  turn,  generally  agree  with  experimental  observations. 

For  radiation  transport  calculations,  however,  results  are  often 
dominated  by  values  of  the  spectral  absorption  coefficient  where  it  is  small  - 
the  Rosseland  mean  free  path  employed  in  radiation  diffusion  theory  is  an 
inverse  mean  of  the  absorption  coefficient  (Ref.  1,  Eq.  36).  Such  values 
of  the  absorption  coefficient  generally  include  many  small  contributions  - 
including  transitions  with  small  occupation  numbers  and  wings  from  distant 
lines.  Experimental  observation  and  theoretical  prediction  of  many  of  these 
contributions  is  difficult.  (Much  of  the  uncertainty  in  the  magnitude  of 
the  N~  cross  section  involves  estimates  of  such  contributions.)  Accordingly, 
much  of  our  effort  this  past  year  has  been  aimed  toward  an  improved  pre¬ 
diction  of  the  absorption  coefficient  where  it  is  small  by  isolating 


individual  contributions  and  attempting  improved  theoretical  estimates 
thereof. 

In  Section  II  we  review  our  definitions  of  radiation  quantities,  and 
Section  III  treats  some  improvements  in  line  width  estimates,  including  a 
discussion  of  the  line  profile  in  the  line  wings.  Section  IV  describes 
our  molecular  emissivity  model,  with  Section  V  summarizing  our  results. 

For  calculation  of  autoionization  cross  sections  ana  for  numerical  con¬ 
figuration-interaction  Hartree-Fock  calculations  the  integration  of  the 
electron-electron  matrix  elements  between  many-particle  states  of  differing 
orbital  occupation  is  a  tedious  exercise  in  Racah  algebra.  Considerable 
simplicity  is  gained  with  techniques  of  second  quantization  (Ref.  6),  our 
treatment  of  which  constitutes  Appendix  A. 

Some  of  this  work  was  included  in  our  semi-annual  report  (Ref.  7)  - 
of  limited  distribution.  We  repeat  it  here  for  completeness  with  our  cur¬ 
rently  "best"  results. 

II.  The  Radiative  Emissivity 

The  total  power  radiated  per  unit  surface  area  by  a  uniform  isothermal 
plasma  of  temperatures  T,  density  p,  may  be  written  approximately  (Ref.  1, 

Eq.  30  et  seq. ) 

J  «  tt  J*  B(uj)  [1  -  e  11  ^  ^3  du>  (vatts/cn2)  (l) 

0 

in  terms  of  the  black  body  function 

B(u))  du  =  ^  UudU  (watts/cm2/sterad.)  (2) 

tt  e  -1 
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Here  the  photon  energy  u>  and  temperature  T  are  in  eV,  u  =  <u/T  (dimensionless), 
and  o  =  1.028  x  lCr  watts/cm  eV  .  L(cm)  is  a  geometric  parameter  character¬ 
izing  the  spatial  extent  of  the  radiating  plasma: 

(i)  Plane  slab,  thickness  a:  L  =  2a(cm) 

(ii)  Circular  cylinder,  radius  r:  L  =  2r( cm) 

and 

(oa)  *  u(u>)  (1  -  e"10/1) 

in  terras  of  the  absorption  coefficient  p(w)  (cm"1)  for  photons  of  energy  u>  in 
air  at  the  temperature  density  P. 

For  a  thick  source,  (p’L  »  1  for  most  oa)  J  =  oT4,  the  black  body 
result  obtains.  It  is  convenient  to  describe  the  energy  twitted  as  a 
fraction  of  the  black  body  emission  at  the  same  temperature.  We  thus  define 
a  dimensionless  emisslvlty  e(T,p,L)  depending  on  tenperature  T,  density  P, 
and  geometric  parameter  L,  according  to 

J(T,P,L)  s  e(0?,P,L)(cT4)  .  (3) 

That  is, 


:(T,P,L)  =  ^  J  <*>3  4s  “^fT  . 


<*) 


(l-e*w/ A) 

If  (n’L)  «  1  for  most  w  the  Planck  mean  opacity  is  a  convenient  result, 

Up  «  J*  tt>3  u(tt>)  du>  (cm"1)  (5) 

rr  T  o 


in  terms  of  which 


The  physics  problem  is  the  determination  of  the  basic  absorption 
coefficient  n(uo)  and  is  described  in  Ref.  1  and  in  suosequent  sections  of 
this  report.  The  numerical  problem  of  the  evaluation  of  the  integral  (Eq..  4) 
over  the  photon  energy  u>  is  the  topic  of  Appendix  B. 

III.  The  Absorption  Line  Profile 

The  results  presented  in  Ref.  1  assumed  the  atomic  lines  to  be  broadened 
only  due  to  the  interactions  of  the  radiating  atoms  with  free  electrons  -  Stark 
broadening.  At  higher  temperatures  or  lower  densities  the  Doppler  width  may 
dominate  the  Stark  width,  and  at  lower  temperatures  and  high  density  the 
resonance  and  Van  der  Waals  atom-atom  interactions  become  impoitant.  The 
results  presented  here  include  both  these  effects,  our  treatment  of  which 
is  described  in  subsections  A  and  B. 

The  Baranger-Griem  classical  path  imr'act  theory  for  line  widths  of 
neutral  radiators  has  been  amply  verified  experimentally  (Ref.  8).  Our 
version  of  this  theory  is  reported  in  Ref.  1.  For  ion  lines,  however,  the 
corresponding  theory  generally  predicts  widths  too  small  by  factors  of  two 
or  more.  Some  recent  theoretical  (Ref.  9)  and  experimental  (Ref.  10)  results 
suggest  the  modifications  described  in  subsection  C,  where  we  include  for 
comparison  our  estimates,  alternative  theoretical  values,  and  experimental  results. 

At  high  temperatures  or  low  densities  the  multitude  of  contributing 
atomic  and  ionic  lines  in  most  all  spectral  intervals  tends  to  preclude  emis¬ 
sion  dominance  at  any  frequency  from  far  wings  of  lines.  At  lower  temperatures 
and  high  densities  the  opposite  may  be  true.  Our  previous  results  for  normal 


density  air  at  kT  =  1.0  eV  provides  such  an  example  (Fig.  1  -  a  more  recent 
version  of  Fig.  8  of  Ref.  l).  The  continuum  absorption  coefficient  for 
photons  of  energy  <i)~  5»0eV  is  =  1.05  x  10  cm"  ;  the  Lorentz  wing  of  the 

9  |i  p  3  Ik 

dominating  NI  line  (2pJ{  rS)  -•  2p  (JP)3s(  P))  at  the  same  energy  contributes 
u(line)  =  2.02  x  10  cm  ,  although  the  line  center  and  width  are  10.33  eV 
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and  3  x  10  eV,  respectively.  That  is,  assuming  a  pure  Lorentz  profile, 
the  continuum  is  dominated  by  a  line  located  over  1000  half -widths  away. 

The  inapplicability  of  the  simple  Baranger-Griem  impact  theory  to 
such  circumstances  is  well  known.  The  possible  magnitude  of  the  error 
resulting  from  ignoring  the  appropriate  validity  criteria  is  a  new  result.  In 
subsection  D  we  describe  come  new  line  wing  approximations  we  intend  to 
incorporate  Into  our  codes  so  as  to  predict  the  absorption  coefficient  in  such 
wings  more  reliably. 

IIIA.,  Doppler  Broadening  and  the  Voigt  Line  Shape 

Due  to  thermal  motions  of  the  radiating  atoms  an  infinitely"  sharp 

spectral  line  at  photon  energy  u>o  takes  on  the  Doppler  shape  (Ref.  11) 

2 

^(u,)  =  exp  [-( - sr-2-)  ]  (7) 

/ff  WD  WD 

where  the  Doppler  width  is  given  by 


M  is  the  mass  of  the  radiating  species. 


•V* 


For  air  atoms: 


v»o  ~ i-2  * io'5  rer 


For  air  molecules:  (diatomic) 

Vwo«  °-85  x  io_5y^v’  * 

The  collision  width  of  Ref.  1  depends  on  the  relative  velocity  of  the 
radiator  and  perturber.  For  electron  perturbers  the  radiator  thermal  veloci¬ 
ties  are  much  smaller  than  electron  velocities  and  the  Lorentz  and  Doppler 
broadening  mechanisms  may  be  independently  treated.  For  atomic  perturbers, 
however,  the  comparability  of  thermal  and  collisional  velocities  should 
couple  the  Lorentz  and  Doppler  broadening  (Ref.  3,  12).  At  the  lower 
temperatures  when  atom-atom  perturbations  are  expected  to  be  important, 
resonance  broadening  is  velocity  independent  (see  Section  B  below),  but 
Van  der  Waals  broadening  varies  as  where  v  is  the  relative  velocity, 

so  a  coupled  approach  may  be  expected  to  modify  the  line  shape  somewhat. 

The  results  reported  here  assume  Lorentz  and  Doppler  broadening 
mechanisms  are  independent,  the  resulting  line  profile  being  the  Voigt 


profile  (Ref.  13). 


V”>  =-37 h 

rr  '  w_  -®  x  +  y  K 


01-0) 


*  =  )  >  »  -  V“c 


■*»  > 


and  Wq  is  the  Lorentz  collision  width  of  Ref.  1.  We  approximate  this 
function  as  follows: 


(i)  (k  <  0.2  or  J  xj  >  5h)  :  (Modified  Lorentz  limit,  Ref. 14) 

y*>)  *  mr [ - 1 — -—5  ♦ - — - -7]  (9») 

^  Wc  1  +  (x  +  kA/5)  1  +  (x  -  h/v£)2 

(ii)  ( |  xj  s  5k  and  k  *>  5.0) ;  (Modified  Doppler  limit) 

<u  -  ok  ^ 

.  / _ °)  2 

bv(»)  -  te  '  ^  [1  +  4(1  -  2  4)]  +  i  f(-UL)}  (9b) 

JZ  WD  h  k  K  K 

where  f(y)  is  the  slowly  varying  function 
tabulated  in  Table  I. 

(iii)  (  |  x  |  <:  5k  and  0.2  ini  5*0)  Tabulated  values  from  Ref.  14 
(reproduced  in  Ref.  1%  pg.  B-l4) 


As  discussed  in  Ref.  1  estimates  of  the  eoissivity  contributions  of 
individual  lines  are  useful  for  identifying  critical  transitions  as  well  as 
providing  a  monitor  to  the  more  detailed  (and  presumably  more  accurate) 
numerical  Integration.  For  a  purely  Lorentz  lire  shape,  the  Ladenburg- 
Reiche  approximation  applies  (Ref.  1,  Eq.  39 )•  For  the  Voigt  profile 
no  such  simple  approximation  obtains.  Recently,,  extensive  numerical 
tabulations  of  "equivalent  widths”  for  Voigt  profiles  have  been  published 
(Ref.  16)  which  should  be  useful  for  detailed  work.  For  our  purpose  a 
simpler  approximation  is  sufficient.  We  employ  the  Curtis  approximation  - 
quoted  in  Ref.  IT  (p.  72)  -  according  to  which  the  eraissivity  of  an  isolated 
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spectral  line  with  Voigt  shape  may  he  approximated 


CLRCD  /,  *D\ 

€V  =  ®LU  +  “  (-  s  =  CLR  1  ”  T  *1 

(lAph)  PpL 


(10) 


€LR  *s  ^adenburg-Reiche  result  and  is  the  isolated  line  emj.ssivity 
for  a  pure  Doppler  profile, 


,  -  ”  <-V 

eD  “  ^  E  - 


n 


n--o  Jn  +  1  (n  +  l)f 


X. 


where 


-  5  <*»  <  *> 


«  4 


w„ 


=  "(iff)  ^L(<»>c) 


flHxT-  (11) 

X^J  J  —  (*D  26^ 


in  terms  of  the  absorption  coefficient  at  line  center,  assuming  purely 

Lorentz  collision  broadening. 

Our  computer  program  employs  the  Curtis  approximation  (Eq.  10)  with 
the  Doppler  approximations  (Eq.  ll).  Comparison  with  our  detailed  line 
integration  results,  in  cases  where  one  expects  an  isolated  line  approximation 
to  be  valid,  indicates  these  approximations  are  generally  valid  to  better  than 
ten  per  cent  accuracy. 


IIIB.  Resonance  and  Van  der  Waals  Broadening 

At  high  density  and  low  temperatures  the  number  density  of  free  electrons 
becomes  sufficiently  small  that  atom-atom  broadening  becomes  important.  The 
theory  of  such  broadening  has  been  worked  out  by  Griem  (Ref.  3) .in  a  perturb¬ 
ation  theory  approach  with  a  "strong  collision"  cutoff  chosen  to  satisfy 
unitarity  (analogous  to  the  Stark  broadening  treatment  of  Ref.  l).  Here  we 
merely  reproduce  his  results  in  the  form  employed  in  our  calculations. 


The  line  absorption  coefficient  for  ground  state  atoms  is  broadened 
by  the  resonance  exchange  collisions  of  the  excited  (final  state)  atom  with 
other  ground  state  atoms.  If  H  is  the  number  density  of  the  ground  state 

D 

atoms,  fR  and  the  absorption  oscillator  strength  and  energy,  respectively, 
of  the  resonance  line  under  consideration,  the  line  width  is 


(g^  and  gg  are  statistical  weights  of  the  lower  and  upper  states,  respectively.) 

Except  for  resonance  exchange  collisions  atom-atom  perturbations  are 
dominated  by  Van  der  Vaals  interactions.  If  J*n  is  the  total  number  density 
of  neutral  perturbers,  Griem  (Ref.  3)  estimates  the  broadening  as  due  to 
strong  collisions  done, 

WV  =  h<v>  Kn  ’ 


where  the  polarizability  entering  PMin  is  estimated  by  a  closure  approximation 
similar  to  Ref.  1. 


_  /Snh5?  x 

Pmn '  (IS^! 


1/5 


-  (f 1/5  [%  «  (M)2  r2]'75. 


^  " 
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Here  r  is  the  mean  squared  atomic  radius  of  the  upper  state  in  atomic  units 
(Ref.  1,  Eq.  22)  and  Ep  is  a  representative  excitation  energy  -  taken  to  be 
the  energy  of  the  strongest  resonance  line.  For  neutral  oxygen  and  nitrogen 
Ep  -  10  eV  =  0.7t  Ryd  agrees  within  20$  with  the  experimental  values,  which 
is  likely  within  the  validity  limits  of  the  closure  approximation  itself. 
With  these  approximations  we  obtain 

vv  -  6.87  x  10-22  s  (-Jf)3/l°C(^)2  r5]  /5  (ev> 

v  11  Me  EP 


^  8.0  x  10' 


■  26  3/l0  *~6‘  2/5 

25  N  (cm  3)(T  „)  (r)  (eV) 


for  air.  Observing  that  7  ~  5v  /2  in  terms  of  the  "effectivu 

number'1  v  of  the  upper  state  (Ref.  1)  and  writing  Nn  =  5*37  x  10^(p/p  ), 

Van  der  Waals;  Wy  s*  6.4  x  Kf  5(p/pQ)  T3/10  v8/5  (eV)  (l4a) 

For  comparison,  we  list  the  corresponding  expressions  for 


Doppler: 


«D  =  *  10  %v  <eV' 


(l4b) 


Resonance:  WRgs  M  5 .6  x  lcf  2(lFgp/p0)  /fl~  (^/(^(eV))  (eV)  (14c) 

4  82 

Adiabatic  Stark;  Wg  »  1.0  x  10“3(Fip/po)T^6  v8/3  (eV)  (l4d) 

where  the  approximations  of  eq.  l4a  nave  been  applied  to  eq.  21  of  Ref.  1. 

F  is  the  fraction  of  the  atoms  in  the  ground  state,  and  F.  is  the  fraction?!. 

£  i 

ionization  of  the  plasma  (assuming  singly- charged  ions  only) .  (All  the 
numerical  values  above  are  for  air  atoms  only.) 


The  results  presented  in  this  report  assume  the  collision  wid:h,  W^, 
to  "be  given  by 

*C  -  Ms  +  V+  V 

where  Wg  is  from  Ref.  1,  (Eq.  13,  above),  and  WR,  where  applicable,  from 
Eq.  12  (above).  The  resulting  Wj,  is  combined  with  the  Doppler  width  as 
described  in  Section  A  above. 

IIIC.  Electron  Impact  Broadening  of  Ion  Lines 

Numerous  attenqxts  have  been  made  to  explain  the  discrepancy  between 
classical  path  impact  theoxy  and  experiment  for  ion  line  widths  (see,  for 
example,  Ref.  18-20).  Although  the  resulting  discrepancies  were  sometimes 
decreased,  they  were  sometimes  increased.  Overall  differences  as  much  as 
factors  of  two  generally  remained. 

In  view  of  these  uncertainties  in  theoretical  prediction  we  chose  to 
follow  the  semi-empirical  approach  familiar  in  astrpphysical  applications 
(Ref.  1  and  Ref.  21,  for  example),  an  approach  Independently  suggested  by 
Griem  (Ref.  22).  The  line  width  is  written 

W1  -  I  <  neT  °)  M  >Av  <»> 

as  a  thermal  average  over  the  electrons  (cm" ^)  of  the  total  inelastic 
scattering  cross  section,  which  we  write 

a,(v)  =  £  CwJv) 


■  ?  f Jk  e)V 

K  *  >  Jk 
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e  is  the  energy  of  the  incident  electron,  uu^  and  f ^  the  energy  difference  and 
f -number,  respectively,  between  states  j  and  k,  and  the  "effective  Gaunt  factor" 
g(co„  ,  e)  is  chosen  to  reproduce  collisional  excitation  experiments.  (Note:  if 
u)..  >  0,  f is  the  absorption  f-numberj  if  io.v  <  0,  it  is  the  emission  f-mtaiber. ) 

jK  JK  JK 

Performing  the  thermal  average  indicated  in  Eq.  15,  we  obtain 


w^(eV)  =  5.22  x  10' 


ne(cm"“) 

/T(eV) 


in  terms  of  the  thermal-averaged  Gaunt  factor 

00 

G  (u>,T)  =  J  e’ydy  g(u>,yT)  . 

0 

Our  previous  results  followed  van  Regemorter  (Ref.  23),  whose  Gaunt  factor  we 
approximated  (Ref.  24,  Eq.  4-142) 


c  0,20 

•=  0.052  In  x 

*  gx  s  0.30  (In  x  -  1.2) 

*  gC(\>tlf) 


(e  <  So) 

(€  >»«,  0  <  In  x<  1.45)  (l8a) 

1  Interpolated  for  (1.45  < 

I  In  x  <  6) 

/  In  x  >  6 


where 


k-k’| 


_  (  c-  )a  ”2 
K  “  ' Ryd'a0 


<  <*’>2  - 


gC(T)t,Tlf)  is  the  Coulomb  Gaunt  factor  (Ref.  25) 

We  interpolate  linearly  vin  In  x)  between  g.^  and  gC  for  (1.45  <  In 

x  <  6). 


Recently,  Bely  and  Griem  (Ref.  9)  employed  t-matrices  for  Mg+  evaluated  by 


close-coupling  calculations  to  evaluate  the  line  width  of  the  resonance  trans¬ 
ition.  They  express  their  results  in  terms  of  an  "effective  Gaunt  factor" 

(through  Eq,  15,  1 6  above)  and  find  results  comparable  to  those  of  van  Regemorter 
near  threshold  and  approximately  a  factor  of  two  greater  than  his  at  higher 
energies.  We  fit  their  results  as  follows? 

g^io,  e)  =  0.160  (®  <*) 

=  gQ  +  ^55  (°*41  -  go)  (0  <  In  x  <  1.65)  (18b) 

s  gj  =  0.13  +  0.17  In  x  (1.65  <  In  x  <  3.0) 

go  =  gBG(u)'u,)  =0.160+1  ^  . 

For  u>  «  c  general  considerations  lead  one  to  expect  g(o>,  e)  “*  gC(Tli,T)f ),  the 
Coulomb  Gaunt  factor.  Thus  a  suitable  extrapolation  procedure  between  g^,  and 
gc  is  necessary.  For  our  calculations,  we  choose 

gpfaj®)  =  ££,(<»>*)  (In  *  <  I.65) 

*(^%3^)sl  +  (~f>^---^)gC(T|1,Tlf)  (1.65  <  In  x  <  3)  (18c) 

=  gC(Tl1^f)  (In  x  >  3) 

In  Figure  2  are  the  thermal -averaged  Gaunt  factors  obtained  from  these  alternative 
approximations,  for  temperature  M?  =  2.0  eV.  Using  these  Gaunt  factors  we  have 
evaluated  the  line  widths  of  some  recently  measured  Nil  lines  (Ref.  10).  The 
results  are  in  Table  II,  including  our  previous  results  and  some  earlier  theor¬ 
etical  values  as  well  as  the  observed  experimental  values.  Comparing  the 


Bely-Griem  (Eq.  l8b)  and  Bely-Griem-Coulomb  (l8e)  Gaunt  factors  of  Figure  2,  it 
is  evident  that  good  agreement  depends  as  much  on  a  suitably  chosen  interpol¬ 
ation  to  the  coulomb  limit  as  on  threshold  details, 

The  Gaunt  factor  of  Eq.  18c  has  been  included  in  our  present  code. 

Figure  3  presents  the  associated  thermal  averaged  Gaunt  factor  at  several 
temperatures . 

HID-  Line  Wing  Absorption 

As  described  in  the  introduction  to  this  chapter,  photon  absorption  in 
the  far  wings  of  strong,  isolated  lines  may  contribute  substantially  to  the 
valleys  in  the  total  absorption  coefficient,  especially  at  the  lower  tempera¬ 
tures  and  higher  densities.  Baranger  and  Stewart  (Ref.  27)  early  considered 
the  absorption  in  windows  between  spectral lines  of  ions  and  introduced  the 
"one-electron"  approximation.  Griem  (Ref.  28-30)  has  extensively  studied  the 
line  wings  of  hydrogen  in  the  one-electron  approximation.  When  the  atom-perturber 
interactions  -  as  measured  by  the  impact  theory  widths  -  are  of  the  order  of  the 
splitting  of  the  atomic  levels  with  the  same  principal  quantum  number  Griem’ s 
hydrogen  results  should  apply.  We  are  principally  concerned  here,  therefore, 
with  isolated  lines. 

Near  the  line  center  the  impact  approximation  is  generally  valid  for 
electrons,  and  the  line  profile  obtains  from  averaging  over  the  large  number 
of  contributing  weak  interactions.  Assumptions  regarding  lack  of  correlation 
among  the  perturbers  during  the  perturbation  "time,"  or  equivalently,  non¬ 
overlapping  (in  time)  of  the  atom-single  perturber  interactions,  allows  the 
average  over  the  product  of  the  perturber  interactions  to  be  expressed  as 
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the  product  of  the  averages  and  the  usual  Lorentz  profile  obtains. 

If  Axi  is  the  displacement  of  the  photon  energy  of  interest,  tu,  from 
the  energy  of  the  unperturbed  line,  j 

Am  =  iv  -  «uij  su>  ■  («i  *ej)  >  v1?) 


and  Wj  is  the  inpact  theory  line  width,  the  Lorentz  wing  profile  is 


b  (to)  — 

J  A  u)  »  Vj  tt(A  w) 


2 


(20) 


where  the  collision  part  of  the  width  is  given  by  (Ref.  1,  Eq.  18) 


f(v)  is  the  Max-, fell  velocity  distribution  of  the  n  (co“^)  perturb  ers  of  mass 

_  C 

m,  X  -  ~  *  j  <k]  r  |  j  >  |  ,  averaged  and  summed  over  magnetic  substates, 

and  A(Z)  is  expressed  in  terms  of  Bessel  functions  of  imaginary  argument 

A(Z)  .  Z2  [K02(|  Z  K^Z)]  .  (32a) 

Comparing  with  Eq.  16,  noting  R.v  =  3(~^“)f we  see  the  function 

JK  n 

a(Z)=J^-A(Z«)  (22b) 

Z  Z* 

plays  the  role  of  an  ’’effective  Gaunt  factor”  for  electron-atom  collisions, 

2 

This  classical  path  impact  theory  is  valid  if  p^in  »  X.  ,  if  kT  »  oUjk  contri¬ 
buting  importantly  to  Eq.  21,  and  A  u>  «  =  (h  x  Plasma  frequency). 

Far  from  the  ?  2  center,  the  perturbation  "time”  t  =  h ft,  w  will  v'Come 

short  enough  that  perturber  motion  may  be  neglected.  Then  the  quasi-static 
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approximation  is  valid.  The  line  profile  may  then  be  described  in  terms  of 
the  static  field  distribution  func -ion  W(F).  For  isolated  states  non¬ 
degenerate  in  orbital  angular  momentum,  the  quadratic  Stark  shift  applies, 
according  to  which  energy  levels,  ej(F),  in  a  field  F  are  shifted  from 
their  unperturbed  values  e  j  by 


* 


where 


°J  '  V  l  1 JX  •  (23) 

The  line  profile  is  then  given  by 

CD 

bj?(»)  «  f  dFW(F)  5 (A  o>  +  C.F2)  (24) 

0  J 

if  only  the  upper  state  is  significantly  perturbed.  (Inclusion  of  the  lower 

state  replaces  by  (C^  -  C^).)  For  the  simple  Holtsmark  distribution  (Ref.  26) 

F  3/2 

W(F)  dF  *  |  (j?)  ~ 


where 


Holts */  \  3  0\ 

'Ul'  4A  u)  'A  (D  ' 


2  ,  4/3  C, 

Au>0  =  (13.62)0*8^)  (-^)  (Ryd) 


for  N  quasi-static  perturbers  (cnf^)  of  mean  charge  2.  The  quasi-static 


distribution  is  expected  to  apply  for  frequencies  (Ref.  29,  Eq.  8) 


(26) 


In  the  intermediate  region, 


the  one-electron 


approximation  is  valid  if  the  impact  approximation  is  valid  at  the  line 
center.  As  discussed  by  Baranger  (Ref.  ll),  the  one-electron  theory  treats 
the  photon  absorption  as  a  three-body  collision, 

e  +  hv  +  (Atom)  e*  +  (Atom)  , 


where  the  spectator  electron  provides  overall  energy  conservation.  Thus,  in 
the  red  wing,  A  id  <  0,  the  electron  must  provide  energy  A  u>,  while  it  absorbs 
A  u)  in  the  blue  wing.  This  results  in  a  line  asymmetry  as  A  ib  •*  H  due  to 
the  population  depletion  of  electrons  with  energy  e  >  A  <t>,  hence  a  suppression 
of  the  red  wing. 

Clearly  such  an  energy-conserving  picture  requires  that  the  interaction 
be  completed  in  the  time  t  ~  h/A  uk  In  terms  of  classical  path  impact  para¬ 
meters,  only  electrons  with  velocities  v  and  inpact  parameters  p  satisfying 
p  <  (~j)  may  contribute,  the  so-called  Lewis  cutoff  (Ref.  31)  •  Our  intention 
in  this  section  is  to  provide  a  formulation  of  the  line  wing  absorption  profile 
which  varies  continuously  from  the  impact  to  quasi-static  limits,  delineates 
the  appropriate  limiting  conditions,  and  incorporates  the  above-mentioned 
line  asymmetry  due  to  electron  population  depletion. 

From  standard  time-dependent  perturbation  theory  one  easily  obtains  the 
line  profile  function  (see,  for  example,  Ref.  3,  Eq.  4-21,  or  Ref.  11,  Eq.  9-23). 
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«  -  i  (lUHU  a)s 

=  =  Rs  £  P„  <  a  |d 3  i  P  >  <  P*  Id  |  a  >  r  dse  “P  {<P  I  U(s,0)  |  p-  >  } 

of 80*  a  "  °  Av 


where  we  have  made  the  usual  classical  path  approximation  and  have  assumeo.  only 
the  upper  state  to  be  significantly  perturbed.  PQ  is  the  density  matrix  -  asvuned 
diagonal  -  for  (unperturbed)  lower  atomic  state  a,  d  is  the  electron  dipole  moment 
operator,  U(s,o)  is  the  interaction  representation  time-evolution  operator  of 
the  atomic  states  p  -  due  solely  to  perturber  interactions  -  and  the  average 
{ . . .}Av  is  taken  over  all  perturber  classical  paths.  Expressing  the  time- 
evolution  operator  in  terms  of  the  time-dependent  potential  due  to  the  classical 
perturber  trajectories,  and  retaining  only  the  second-order  (lowest  non-trivial) 
term,  one  obtains  (after  interchanging  orders  of  integration) 

go)  -  ff  Be  E  pa<a  |  4  |  p  >  <p-  |  a,  |  o£(4)  jate  T  0,9  JttV*  1 

<*&•<>  f»2  ■  o  0 


Jt2 


{<  p  i  v(t, +tg) !  p2  >  <  pp  1 v^)  i  p*  >  }AV  . 


In  the  one-electron  approximation  the  average  is  to  be  taken  over  all  classical 
paths  of  a  single  perturber.  Assuming  the  statistical  distribution  of  perturbers 
is  stationary 

{  f(t  +t)  f(t)  }Ay  =  {  f(r)  f(0)  }Ay  , 
we  obtain  (Ref.  29,  Eq.  12) 

1  Be  T  »  <0<)  aa  1  e>  <P'  I  1  )\ 


<*p/v — ap* 


(  <  p !  v(tx)  1  p2  >  <  p2 !  v(o)  |  p»  >  }Av 


and  'll  i  lil  Qr*  m 
—  -  r  xap  -ap»  * 


Substituting  the  classical  potentials 

2  3  i  , 

<  p  |  V(t)  i  p‘  >  =  e2£  O  I  Kv  I  3’  > 

v=l  I  r(t)  I 


'r(t)  =  P  +  v(t-t0),  p.v  =  0 


and  performing  the  one-electron  average 

00  d2fi  pMax  * 

ne  J*dv  f(v)  |  2npdpj  vdtQ  (., 


one  easily  obtains 


f  <  3  I  v(t1)  1  32  >  <  p2  j  V(0)  |  p‘  >  } 


n  e4  ® 
e  (» 


Max  ® 


[e  <p1rvIp2><p2|rvIp»  >]. 


Av  =  -J—  Jf(v)dv  Jjnpdp  Jvdtc 
o  pMin  -» 

[p2-v2t0(trt0)] 

T1 - rwriT^ 

[p2+v"t02] 


Then,  with  substitutions  ^ 

ZaP2  =  ^v  ^aPg"^^  the  time-integrations  result  in 


L(**»)  =  J  Re  ]C  P, 


<  a  |  d  |  p  >  <p*  j  a.  j  oc  >  n  e4 


app’a  a  *  3?v 


T~  A?  <  P  I  EJ  ^"  >  <  3”  I  Rj  3’  > 


w  raeuc  _ 

£  f(v)vdv  J2npdp  (-g-g  ■)  C  [^L  “>)] 


where  C(X)  =  A(X)  +  iB(X)  is  defined  in  Ref,  1.  (Note,  i.(x)  is  even  in  X.) 


In  the  line  wings,  we  approximate  6*  =  0  and  extract  the  line  shape 
function 

oo  pMax 

b(“)  *  7177 1  "e  J'' av  Hv)  hT  >S  k  V  J"  f*  A(-&  («■”»>)  (S7) 

W  °  pMin 

which  compares  exactly  with  Eq.  20,  21  in  terms  of  the  same  A-function  (Eq.  22a) 
but  evaluated  with  a  frequency-dependent  argument, 

Mik  ■  “^u^dk  - A  “  -  v.  - 

Equation  (27)  is  our  basic  result;  it  remains  to  evaluate  the  p-limits. 

The  completed  collision  criteria  requires  that 


Max 


=  p  =  Min  { -Si  ,  J  s 
2  1  (u  *  Auu  } 

P 


hv 


(28) 


The  unitary  condition  requires 


o  p  p  ,  ^PMin  ^Min  , 

PMin  S  3  *  \  V  1  A  ClvT  (A  +  i  X-ftr  (A  I 

From  the  definition  of  PMax>  the  largest  value 

PMin  A  “"<•>  ft 

(A  uj-o) j^)  =  ( — - — *=)  <  1,  and  from  tabulated  values, 
|  A(x)  +  i  B(x)  |  —  1  for  0  <  x  <  1,  so  we  take 


4n5 


(29) 


This  choice  of  p„.  evidently  satisfies  the  validity  condition  for  the  classical 
Min 

path  approximation. 

Near  the  line  center,  A  uj  "  0,  pg  -*  »  and  the  impact  theory  approximation 
obtains.  Thus  the  validity  conditions  for  impact  theory  must  include  Am  «  u>^, 
where  k  is  the  state  lying  nearest  in  energy  to  the  state  j  and  connected  with 
it  by  a  dipole-allowed  transition.  Collecting  these  results,  our  frequency- 
dependent  "line  width"  becomes 
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1 


WjW-ij  ^5  ne  |  T  f(v)  l  Ejk  /  A  [  (A  “  '  V)] 

mo  K  "l 

"  %  a  ( Vo3)  w  (Ryd) 


(30) 


where 


W(Ryd)  -E  R.J  dye  y  J  ~  A(ZX) 
k  v  Z(y)ZX 

yM 


z2  •  ^ — 2ui — ■/  ’  Zl  ■  y  Z2  ’  yM  ■ 


/hAjjx 


Perturbers  of  energy  e  <  yM  MP  have  been  removed  from  the  one-electron  theory.  Ac¬ 
cording  to  Eq.  26,  they  should  be  added  to .the  quasi  static  profile. 

Generalizing  the  effective  thermal  Gaunt  factor  (Bef.  1,  Eq.  21  and 

Fig.  1) 


CO 


a  (X,  V)  *  ;  dte_t  a(x/t)  -e~YA  a  (±) 
YX 

we  write 

/  T 

W(Ryd) 


(31) 


(32) 


The  function  a(X,  y)  is  presented  in  Fig.54  for  0  <  y  <  1.8.  Near  the 
line  center.  Aid  -*  0  and  a(X,0)  reproduces  0t(X)  of  Ref.  1.  Far  into  the  line 
wings,  y  -*  1. 

To  incorporate  the  red  wing  electron  depletion  effect  into  our  profile, 
the  classical  path  theory  is  no  longer  adequate,  as  it  ignores  the  change  in 
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perturber  energy  during  the  interaction.  It  is  well  known  from  coulomb  excit¬ 
ation  theory  (Kef.  32)  that  a  classical  coulomb  Gaunt  factor  with  sysmetrized 
energy  argument  e  =  ^  (e^  +  e^)  frequently  provides  an  adequate  representation 
of  the  quantum  mechanical  result.  We  similarly  symnetrize  our  Gaunt  factor, 
obtaining 


The  strong  collision  impact  theory  term  must  be  similarly  treated.  We  find 

2/3  -yM  2/3 


w®(u»)  »  2nd  (i 


W’fW  [i  V  (a%)]  r  <!>  (34) 


and  the  total  one-electron  width  is  the  sum  of  Eq.  3^  and  30,  with  W  from  Eq.  33. 
Finally,  the  number  density,  N,  in  the  quasi  static  term  (Eq.  25)  must  be  in¬ 
creased  by  the  number  of  electrons  removed  from  the  one-electron  term 


N 


N  +  n 


J 

0 


'M 


2 

7sr 


y1/2  e“y  dy 


(35) 


With  these  expressions  we  have  obtained  a  computationally  feasible  formul¬ 
ation  of  the  absorption  line  profile  varying  continuously  from  the  impact  theory 
line  center  to  the  quasi-static  far  line  wing.  These  equations  have  not  yet 
been  incorporated  into  our  radiative  mean  computer  programs,  but  will  be 
shortly. 
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IV, 


The  Molecular  Absorption  Coefficient  and  Emisslvity 
For  temperatures  below  1ST  =  1.0  eV,  absorption  by  molecular  species 
becomes  of  primary  importance.  Avilova  et  al,  (Ref.  4)  recently  siuanarized 
molecular  contributions  to  the  air  absorption  coefficient.  Comparison  with 
earlier  vork  (Ref.  51  and  53)  Indicated  some  large  differences,  toward  the 
resolution  of  which  most  of  the  work  here  described  has  been  directed.  In 
almost  all  cases  our  results  will  be  found  to  oe  in  general  agreement  with 
Avilova  /Ref.  4). 

Ir  subsection  A,  below,  we  define  our  notation  and  write  down  our 
representation  of  the  absorption  coefficient  of  a  single  molecular  line.  We 
have  employed  both  the  "smeared  band"  model  of  Avilova  (Ref.  4)  and  a  more 
detailed  three-branch  "PQR -model,"  as  described  in  subsection  B.  The  molecular 
bands  we  have  studied  and  the  associated  data  constitute  subsection  C,  and  in 
subsection  D  we  describe  some  emission  band  models  we  have  studied.  Finally, 
in  subsection  E,  we  briefly  describe  our  molecular  absorption  coefficient 
computer  program  ABSMOL. 

Our  notation  closely  follows  Herzberg  (Ref.  38)  and  Nicholls  (Ref.  52), 
and  we  shall  frequently  follow  the  convention  of  choosing  wave  number  k  as  our 

4 

energy  variable,  noting  that  a  wave  number  10  cm  corresponds  to  an  energy  1.240  eV. 

IVA.  The  Absorption  Coefficient  of  a  Molecular  Line 

Molecular  states  are  labeled  according  to  their  electronic  state  e,  the 
component  of  electronic  angular  momentum  along  the  internuclear  axis  Ay  and  the 
rotation  and  vibration  quantum  numbers  J  and  v,  respectively.  Upper  states 
ere  singly  primed  (')  and  lower  states  are  doubly  primed  ("). 
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We  consider  a  transition  from  a  lower  i.tate  i  £  (v",j,,,e",A") ,  one 
of  states  of  energy  e^,  to  an  upper  state  j  £  (.v’jj^e', A')  of  energy 
en  and  degeneracy  dfi.  The  f-number  for  such  a  transition  is  (Ref.  38  ,  52) 

e_ 


J’A’  Re  ^v'v"^ 


f  _  1  ,*mtu  J’A’  e 

m-n  “  3  'Ryci/  Vv"  bJ"A"  d 

m 


(36) 


in  terms  of  which  we  may  express  the  absorption  coefficient  as 

>Wk)  =  nr0Nmfm-nbmn(H)  * 

2,  2 


(37) 


(r  =  e^/mc11  =  2.8l8  x  lCf1^  cm) 


N  (cm  ^  the  number  density  of  molecules  in  the  (lower)  state  m,  and  b  (k)  (cm) 


mu' 


is  a  line  shape  factor,  normalized  according  to 

J  **  W**  =  1  ‘ 

o 

In  eq.  36,  =  (en  -  em)  is  the  transition  photon  energy,  „  the  Franck- 

Condon  factor,  Sj„^„  the  Honl-London  factor,  and  Re2(rv,v„)  the  electronic 

transition  moment  in  atomic  units.  In  writing  down  eq.  36  we  have  assumed 

the  validity  of  the  Born-Oppenheimer  approximation,  and  we  shall  further 
2 

assume  Rg  to  be  a  smoothly  varying  function  of  its  argument  r^,^,,,  the 


r-centroid  (Ref.  52) 


To  relate  Nm  to  the  species  number  density  N,  we  introduce  the 

partition  function  Q  and  decompose  the  state  energy 

into  electronic,  vibrational,  and  rotational  parts, 

•n  =  *e"  +  V  *  *j"  . 


or,  in  wave  number  units, 

1 


*_  =  T „  +  G  ,i  +  F_it  (cm  )  . 


he  m  "e"  _v"  «r 

T  „  is  the  electronic  energy  of  the  state  (measured  to  the  minimum  of  the 
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potential  curve)  and  we  approximate  (Ref.  38) 


Gv  s  G(°)  +  G0(v)  =  G(o)  +  wQv  -  aJ0xovt  +  (uoyQv3  +  . . . 


Pj  =  BVJ(J  +  1) 


Bv  *  Be  "  ae(v  +  & 


d  -*_/kT 
Nffi  -  N(-|)  e  * 


dm  =  (2J»  +  1)  de„  ,  de#  =  (2SM  +  1)(2  -  «A„^0) 


and  the  usual  approximations  -  assuming  kT  »  hcBy  -  result  in 


Q  - 

%  -  9/By  . 

where  Q  =  kr/hc  is  the  temperature  in  cm"1.  Gilmore  (Ref.  40)  tabulates 

„  ,  .  de"  “«e"/kT 
fe„«  - 

for  most  molecules  of  interest  to  us.  Writing  Ne„  *  fe„(T)N,  we  finally 


obtain 


B  ..  -u>"/e 

Nm  =  V  (“f)(l-e  °  )  (2j"  +  1)  e 


"[G  (v")  +  Ft„]/8 


Combining  the  above  expressions  our  final  form  of  the  absorption 


coefficient  is 


,  ,  "roN  Bv’Vi  -CGo(v)+FJtt]/e 

iwW  T  Mt*^1  "  e  '  * 


-  A wN  .  oJ'A1  fRe2(rv*v”\  v  ,  \  ,_-lx 

x  W  Vv"  Sj»A”  C - dT - J  V(h)  (ca  5 


-25- 


[Note;  Many  authors  (ex.  Ref.  4) include  the  electronic  statistical  factor 
de„  in  their  definition  of  the  electronic  transition  moment,  so 
care  should  be  t£};.en  when  comparing  numerical  values.] 

We  write  the  transition  energy  variously  as 

Kmn  =  ^e1  “  Te"^  +  ^Gv‘  “  Gv"^  +  ^FJ‘  ”  FJ'^ 

-  V,  +  -  °c(v"”  +  (FJ-  -  V 

=  Kylyll  +  (^"jt  “  Fjll) 

where  koq  is  the  energy  of  the  transition  [e^v1  -  o,j'  =  o,A'] 

«— •  [e",v"  =  o,J"  =  o,A"],  the  so-called  "o-c."  transition. 


IVB.  Molecular  Band  Models 

In  applications  we  are  generally  interested  in  an  averaged  molecular 
absorption  coefficient  -  averaged  over  an  energy  interval  A  h  much  greater 
than  a  molecular  line  width  W  and  generally  large  enough  to  contain  many 
individual  lines.  That  is,  for  a  particular  pair  of  electronic  states  (e* ,e") 
we  define 


e'e 


..(h)  = 


JL 

Ah 


h+Ah 

J  d h' 


v"J" 

v’J' 


=  £ 
v"J" 
v’J' 


r  mn  n 

Lb  Ahj 
mn 


(40) 


summed  over  an.  values  of  v’VjV'jJ’  for  which  is  in  the  interval 
Ah  about  k.  To  this  end  the  (J*  -  J”)  selection  rules  must  be  considered. 
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In  general,  J*  =  J",  J"  +  1, leading  to  the  three  "branches  (Ref.  38,  Eq.,  IV, 
32-34): 


"Q-branch"  (J’  =  J”) 
l,K-branch"  (j'  =  J"  +  1) 
MP-branchn  (j»  ■  J”  -  1) 


Q  / 

"tan  -  kv,v„  +  (Bv,-Bv,^ 

Kmn  =  Hmn  +  2Bv' 

4  =  4-  2v  (J"+1> 


(«) 


For  electronic  angular  momentum  A  ^  0,  details  of  the  A  -  J  coupling 
may  lead  to  further  splitting  of  these  branches.  For  simplicity  we  have  ignored 

ft 

such  coupling;  our  Honl-London  factors  are  those  of  a  rotating  top  and  are  tabul¬ 
ated  by  Herzberg  (Ref.  38,  p.  208).  Care  should  be  taken  with  Herzberg' s  formulae , 
as,  for  A  A  f  0,  they  must  be  multiplied  by  a  factor  of  two  in  order  to  satisfy 
the  sum  rule  (Ref.  52) 


J"  A” 

IS  -  (2J“  +  1)  . 

J»  A* 


Explicitly  writing  the  rotational  sum  terms  of  Eqs.  39  and  40: 


v’v" 


r  Hmn  qJ’A* 
j, T„  Ah  SJ"A" 


-Bv„j"(j”+i)/e 

e 


1 42) 


(43a) 


The  maximum  value  of  J"  we  choose  to  be  the  smallest  of  the  following: 

Ji  *  t  5^  C  -  o-(o)  -  V"  +  w”2  ]J1/2 

(and  the  same  with  (')  replacing  (*')) 

J2  =  5.o  ye/Bv« 
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is  effectively  an 


refers  to  the  Morse-Potential  dissociation  limit,  and  Jg 
occupation  number  cutoff. 

Equations  ( 39-^3)  explicitly  define  our  "PQR-Branch"  model.  All  the 
molecular  absorption  coefficients  presented  in  this  report  were  generated  by 
this  model. 


kT*  *5  « i 

At  higher  temperatures,  9  =  r—  (  ~  lCr  cm  )  becomes  much  greater  than 


Bv„  (  ~  1  cm*1),  and  we  clearly  wish  to  replace  the  sums  over  J'  and  J"  in 
Eq.  43  by  integrals.  Since  J  ~  ~  10  for  T  =  1.0  eV,  we  are  inter¬ 

ested  primarily  in  Honl-London  factors  for  J  »  1  and  J  »  A  (Ref.  38,  p.208). 

R-branch  f-branch  P-branch 


AA  =  0: 

=  +1: 
c  -1; 


S1)  ~  J 


J/2 

j/2 


SJ  ~  j 
J 

J 


J/2 

J/2 


These  approximations  exhibit  the  familiar  effect  (Ref.  38,  pg.  20?) 
that  the  Q-brancli  intensity  i.s  approximately  twice  the  R-  or  P-branch 
intensity  (for  AA  /  0) .  Further,  the  p-  and  R-branch  energies  differ 
from  the  0-branch  energy  only  b ./  a  tem  of  order  Bv„/,r"AB,  where 
AB  -  (By,  -  Bv„) .  For  the  principal  air  bands  this  factor  is  less  than 


10/j",  resulting  in  an  energy  displacement  of  order  ten  percent  for  J"  ~  j. 
For  these  reasons,  plus  computational  convenience,  high  temperature  models 
usually  assume  '^-branch  energy  dependence  only  (Ref.  4)  • 

Assuming  a  Q-branch  energy  dependence  defines  our  "  Smeared  Band”  Model. 


mn 


=  K. 


v’v" 


ABJ"(j"-tl)  j 


/-»0 


J 1  A1 

Using  Eq.  42  £  ®j"A"  '  (2<J"+l)  and  expressing  the  sum 

J 


in  Ec±.  *t3a  as  an  integral,  we  obtain 


v'v"  =  AkTB  I  3x'H 


k-*Ak  *B  ,,(— 


v'v 


AS 


0/6 


for  (K-Jiv,v„)AB  *  0 


(43b) 


-  0  otherwise. 

Our  final  expression  for  the  Smeared  Band  averaged  absorption  coef¬ 


ficient  then  becomes 


-</e 


V 


^(k)  -  “3“  «aQ  Ne”(1  ~  e  )  s  (  e'^'Vv’^d  J 

y  fv  1  g,r 


R  2  -G(v")/Q 


*j/vtv..(  era"1)  (44) 


V  V 


f 4tT  2  ~  ra  1(~— 18  2. 

(— ora  =  2.68  x  10  cm  ) 

3  o 

which  agrees  with  the  result  of  Ref.  4. 

We  note  that  the  sum  rule  guarantees  the  total  oscillator  strength 
conservation;  our  Q-branch  model  errs  only  in  an  energy  redistribution  of 
oscillator  strength,  an  effect  of  less  importance  at  higher  temperatures  and 
larger  averaging  interval  A  h. 

All  the  band  systems  we  have  treated  have  been  evtluated  in  both  the 
"FQR-Branch"  and  "Smeared  Band"  Models.  For  detailed  comparison  with  experiments 
of  high  energy  resolution  only  the  FQR-Branch  model  is  suitable.  Furthermore, 
for  the  purpose  of  developing  energy-dependent  distribution  functions  of  the 
line  strength  -  for  our  band-emissivity  and  opacity  studies  -  the  more  detailed 
line-by-line  treatment  is  necessary.  For  evaluation  of  the  mean  absorption 
coefficient,  however,  averaged  over  an  interval  Ah  >  0.1  eV,  the  results  of 
the  two  models  are  essentially  in  agreement  for  temperatures  T  >  C000°K. 


IVC,  Molecular  Band  Systems  Contributing  to  Air  Radiation 

The  molecular  species  we  have  examined  are  V  V  N2+  and  NO.  In 
Table  II  we  list  each  molecular  band  ve  have  studied,  its  spectroscopic 
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designation  and  name,  and  the  energy  interval  in  which  the  band  radiates.  This 

2 

table  also  includes  our  chosen  values  of  the  electronic  transition  moment  R 

e 

(Eq.  36),  and  the  source  of  our  Frank-Condon  factors  and  r-centroids.  The 
oscillator  strength  is  listed  for  the  Rydberg  states  of  NO  and  finally  the 
energy  averaging  interval  At  (Eq.  40),  determining  the  spectral  details  of 
our  average  absorption  coefficient.  The  maximum  value  of  A  x  used  in  our 
calculations  was  0.10  eV,  chosen  large  enough  to  resolve  the  vibrational 
structure,  generally  of  the  order  0.124  -  0.248  eV. 

The  spectroscopic  constants  T  ,  m  ,  m  x  ,  my  ,  B  ,  a  ,  r  depend  only 

6  6  6  6  •fir  6  6  6  6 

on  the  electronic  state  e  and  the  values  we  used  are  listed  in  Table  IV,  with 
the  source  for  these  constants. 

In  Table  V  we  list  the  Line  Planck-Mean  per  molecule  for  Ng ,  0g,  NO  and 
Ng+(in  units  of  10"1®  cm2)for  the  temperatures  4000,  6000,  8000,  10000,  12000, 
18000,  24000  (°K). 


1. 


a) 

d) 

c) 

a) 

(a) 


Ng  Systems 

For  the  Ng  systems  we  include  the  bands: 
Transition 

A  3eu+-  B3rTg  (First  Positive) 

B  -»  C  \u  (Second  Positive) 

X  V;  -  b‘  ( Birge-Hopf  ield  l) 

o  ^ 

X  h*  -  b  \  (Birge-Hopf ield  2) 

o  ^ 


Ng  (First  Positive) 


Spectral  Inverval  (eV) 
0.4  -  2.7 
2.2  -  4.8 
5.0  -  13.3 
5.9  -  13.1 


For  this  band  system  we  used  the  Frank-Condon  factors  of  Ref.  37*  v'  = 


0-17,  v"  =  0-13.  The  values  of  Rg  were  calculated  from 


R  2  =  0.500  -  2.95 x  1-8.73  x2  (l.l6  A  <  r  <  1.61  A) 
e 

where  X  =  (ry!  yl,  -  r^,,)  in  terms  of  the  r-centroid  rv,  y„  (Ref.  52)  and  re,i, 
the  equilibrium  internuclear  distance  of  the  lover  state.  The  coefficients 

p 

of  the  R  fit  were  obtained  from  Ref.  47  and  the  r-centroid  data  from  Ref.  36. 
e 

The  resultant  cross  sections  agree  well  with  the  results  of  Avilova 
(Ref.  4).  (The  values  given  earlier  (Ref.  53)  are  erroneously  higher.) 


(b)  Ng  (Second  Positive) 

We  used  the  Frank-Condon  factors  of  Ref.  37,  v’  -  0-4,  v"  =  0-12.  lhe 
2 

values  of  R  were  calculated  from 

G 

Re2  =  4.0  -  41.20  x2  +  10.1  x4, 
where  X  =  (ry,Y„  -  re„). 

2 

The  coefficients  in  Re  were  obtained  from  Ref.  49,  and  the  r-centroids  from 
Ref.  35. 

(c)  Np  (Blrgo-Hopfleld  l) 

The  Frank-Condon  factors  are  from  Ref.  37,  v*  =  0-6,  v"  =  0-28.  The 
2 

value  of  R„  were  obtained  from  Ref.  4. 

e 

(d)  No  (Blrge-Hopfleld  2) 

We  used  the  Frank-Condon  factors  of  Ref.  31,  v*  »  0-7,  v"  =  0-28.  The 
2 

value  of  R  came  from  Ref.  4. 
e 
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2. 


Op  System 

In  the  figures  presented  here  (Fig.  11-17)  we  have  included  only  line 
transitions 


Spectral  Interval  (eV) 

X  “*  B  ( Schumann-Runge)  1.0  -  7.0 

o 

Hie  Frank-Condon  factors  were  obtained  from  Nicholls  (Ref.  33)  for 
v*  =  0-21  and  V  =  1-21  and  for  v*  =  0-21,  v"  =  0  we  used  the  values  of  Jarmain 
(Ref.  3*0*  The  r-centroids  and  coefficients  used  in  calculating  Rg  were  obtained 


from  Ref.  45. 

Re2  =  3.0625  -  5.25  X  +  2.25  X2 


where  X  =  rv,v„  -  re„  . 

For  the  energy  interval  >7*0  eV  we  used  the  Schumann-Runge  continuum 
results  of  Avilova  (Ref.  4)  in  our  eadssivity  calculations.  Our  results  for  the 
Schumann-Runge  band  agree  quite  well  with  those  of  Avilova  et  al.  (Ref.  4),  but 
differ  from  those  of  Churchill  (Ref.  51)  whose  data  on  Frank-Condon  factors  was 
incomplete. 


3.  Np+  Systems 

In  this  system  we  Included 

Transition  Spectral  Interval  (eV) 

a)  X  V*  -  B  %+  (First  Negative)  1.5  -  4.9 

b)  X  -»  A  2tt  (Meinel)  0.5  -  4.7 

o  ^ 

(a)  N*  (First  Negative) 

Kxe  Frank-Condon  factors  came  from  Ref,  37,  v*  =  0-20,  v"  =  0-20.  The  value 

2 

of  R  was  obtained  from 
e 
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F.e2  =  0.9  +  22.14  X2  +  136.16  X4 

—  2 
where  X=(r,  „  -  r  !t).  The  coefficients  of  R  were  obtained  from  Ref.  49 
'  v’v  e  e 

and  the  r-centroids  from  Ref.  36. 

(b)  (Meinel) 

¥e  used  the  Frank-Condon  factors  of  Ref.  37,  v*  =  0-20,  v"  =  0-20. 

2 

The  value  of  R  was  obtained  from  Ref.  44. 
e 

4.  NO  Systems 

This  system  was  divided  into  two  separate  groups 

GROUP  A 

Transition  Spectral  Interval  (eV) 


a) 

X  2tt  -  B  2tt  (p) 

1.0  -  7.8 

b) 

X  2n  -  A  V  (y) 

3.r  -  7.9 

c) 

X  2tt-  C  2tt  (6) 

4.0  -  7.8 

a) 

X  2n-  D  V  (3) 

4.0  -  8.5 

GROUP  B 

This 

system  consists  of  the 

Rydberg  series  and  the  Ogsva  2  band  system. 

GROUP  B 

Transition 

Spectral  Interval  (eV) 

a) 

K  2tt  -  D 

0.920  -  1.370 

b) 

H’  2tt  -  C  2rr  (Ry) 

\ 

1.245  -  1.445 

c) 

H  2n  ~  C  2tt  (Ry)  > 

(Feast-Heath-Lagerquist)  1.22  -  1.445 

a) 

F  2A-  C  2tt  (Ry-) 

0.920  -  1.220 

e) 

H1  2n  -  D  (Ry) 

(Feast-Heath)  1.120  -  1.345 
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Transition 


GROUP  B  (Cont*  d) 


Spectral  Interval  (eV) 


f) 

H  %+  -  D  V  (Ry) 

(Feast-Heath) 

1.070  -  1.295 

g) 

D  -  A  2S+  (Ity) 

(Feast  l) 

1.020  -  1.245 

h) 

E  -  C  2tt  (Ry) 

0.920  -  1.295 

*) 

C  2n  -  A  V  (Ry) 

(Heath) 

0.920  -  1.070 

J) 

4  4  - 

a  tt  -  b  £  (Ogawa 

2) 

0.470  -  1.970 

(a)  NO  O) 

The  Frank-Condon  factors  were  obtained  from  Ref.  37,  v*  =  0-19,  v"  = 

■> 

0-23.  We  calculated  Rg  from 

R2  =  0.03844  +  0.124  X  +  0.1  X2 
e 

where  X  =  (rv,v„  -  re„). 

2 

The  coefficients  for  R  were  obtained  from  Ref.  49  and  the  r-centroids 

e 

from  Ref.  36. 

(b)  NO  (y) 

For  this  band  system  we  used  the  Frank-Condon  factors  of  Ref.  37  for 

V*  =  0-5,  v"  =  0-23.  For  v*  =  6-23,  V”  =  0-23  we  used  our  calculated  values 

(subsection  E  below).  We  fit  the  data  of  Callear  et  al.  (Ref.  60): 

Re  *  0.350  (r  <  1.0  A) 

«  -19.90  +  39.0  r  -  18.75  T2  (1.0  <  r  <  1.04  a) 

=  -334.267  +  924.768  r  -  850.2  r2  +  260.0  r3  (1.04  <  r  <  1.14  A) 

=  0.250  (r  >  1.14  A) 

The  resulting  f-numbers  agree  with  recent  experiments  (Ref.  48  (1967))  to  within 
ton  percent. 


-34- 


(e)  NO  (6) 

For  this  band  system  we  used  the  Franck-Condon  factors  of  Ref.  37  >  v*  =  0-4 , 

p 

v"  =  0-23.  The  value  of  Rg  came  from  Avilova  (Ref.  4). 

(d)  NO  (e) 

The  Franck-Condon  factors  were  obtained-  from  Ref.  37  >  v*  -  0-7>  v"  =  0-23. 

2 

The  value  of  Rg  was  obtained  from  Avilova  (Ref.  4),  and  the  resulting  f-numbers 
agree  within  ten  percent  with  experiment  (Ref.  59)* 

Group  B 

(a)  -  (j)  NO-Rydberg  Series 

For  the  NO  (Rydberg)  series  the  Franck-Condon  factors  are  from  Ref.  37  and 
our  calculations.  We  use  Wray’s  experimentally  fit  oscillator  strengths  (Table  4 
of  Ref.  42). 

Our  calculated  values  of  the  radiation  intensity  per  particle  as  a  function 
of  wave  length  for  the  Rydberg  series  agrees  well  with  the  values  presented  by 
Wray  (Figures  9-14  of  Ref.  42). 

( j)  NO  (Ogawa  2) 

For  this  band  system  we  used  the  Frank-Condon  factor  v*  =  0-22,  v”  =  0-23 

of  Nicholls  (43)  extended  by  our  calculations.  The  value  of  R  2  was  obtained  from 

e 

Avilova  (Ref.  4). 

IVD.  Molecular  Emlsslvlty 

The  averaged  absorption  coefficient  of  Eo.  40  above  is  fed  into  the  atomic 
radiation  code  RATRAC  (Ref.  l)  as  an  additional  contribution  to  the  continuum 
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absorption  coefficient.  In  this  approximation  we  neglect  the  increased  self¬ 
absorption  of  the  individual  molecular  lines  in  the  line  centers  where  the 
absorption  coefficient  is  much  greater  than  4.  To  attempt  to  correct  for  this 
omission  we  exploit  the  high  density  of  the  molecular  lines  and  the  approximately 
constant  width  (Ref.  54);  we  treat  the  band  emissivity  statistically  (Ref.  55,  56). 

In  our  band  models  rotational  lines  originating  from  a  given  band  (v*  ,v’:) 
are  distributed  with  a  prescribed  spacing  and  a  strength  exponentially  decreasing 
from  the  band  head.  When  |I(k)  is  large  and  A  k  not  too  small  several  bands 
generally  contribute  in  ary  interval  Ax.  In  most  of  our  numerical  applications 
A  h  =  807  cm-1  (0.1  eV)  and  we  typically  find  500-4000  lines  contributing  in  an 
interval,  t..a  large  number  generally  applying  to  intervals  where  4  is  large  and 
self -absorption  more  important. 

Considering  the  extensive  variation  of  the  Franck-Condon  factors  with 
(v*  ,v")  and  the  large  numbers  of  (v'  ,v")  contributing  to  a  given  Ax,  ve  assume 
the  n  individual  lines  in  Ax  may  be  treated  as  if  uniformly  distributed  in 
energy  and  continuously  distributed  in  strength,  with  the  probability  of  a 
strength  S  in  dS  being  P  (S).  Our  molecular  absorption  coefficient  computer 
program  ABSM0L  (subsection  E.  below)  evaluates  a  histogram  of  the  line 
strength  distribution  about  the  mean  line  strength  in  every 
energy  interval  A  h.  From  such  a  detailed  energy-dependent 
line  strength  distribution  function  we  expect  to  be  able  to  obtain  more  re¬ 
liable  estimates  of  band  emissivities  and  opacities.  This  work  is  still  in 
progress. 

For  the  short  path  lengths  and  lower  densities  of  interest  here,  however, 
the  molecular  line  self -absorption  is  not  a  large  effect,  and  we  have  approxi¬ 
mated  the  line  strengths  as  being  exponentially  distributed.  That  is, 
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P(S)  =  s"1  exp( -S/Sq)  in  terms  of  a  mean  line  strength  Sq(h), 

S0  M  ~  =  4k  ?(k)  (45) 

H 

bbe  Juesn  line  spacing.  Comparison  with  our  detailed  histograms  Indicates 
this  approximation  to  be  generally  satisfactory  for  the  majority  of  the  lines, 
but  a  small  number  of  lines  -  of  order  10  -  20#  of  the  total  -  appear  to  lie  in 
a  distribution  roughly  symmetric  about  a  strength  a  few  times  the  mean  SQ.  As 
these  lines  are  most  strongly  self-absorbed  the  exponential  distribution  can 
only  r  'vide  a  rough  indication  of  the  true  self -absorption  corrections. 

For  a  known  line  strength  distribution,  P(S),  and  a  uniform  distribution 
in  energy,  a  simple  approximation  for  the  emissivity  obtains  (Ref.  56). 

Ac(k)  =  (-^)  J  dH*B(K')[l  -  e’*1  (k  )L] 

O0  H 

_  -W  (L)/d 

.  -  ( — if)  B(k)Ah  [1  -  e  k  h]  (45) 

06 

where  V?k(l)  is  the  average  equivalent  width  of  a  single  line  in  A*, 

Wx(L)  =  J  dS  PH(S)  J  d*»  [1  -  e"Sb(H*)L]  # 


For  a  I^rentz  profile  of  width  W  (Ref,  1,  eq.  11)  and  an  exponential 
strength  distribution,  the  result  is  easily  obtained, 


•  Sb(x')L 
(L)  =  {  ^(h*) 


dn’ 


Sq(k)L 
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The  molecular  emissivity  in  the  band  Ah  is  then  approximately 


Ac(k) 


ttB(k)  ,  f,  r_ _ _ 

— Ah  11  -  exp[-p==™===~- 

ae  ,/i+(  wnt,(K)L/,T 


3} 


(*7) 


Treating  the  molecular  emission  as  an  additional  contribution  u(h)  to  the 
RATRAC  continuum  (Ref.  l)  results  in  an  emissivity 


Ac  (h)  =  — r  B(k)  Ah  (1  -  exp[-n(H)L]}  . 
C  00 


(48) 


It  is  interesting  to  consider  Eq.  47  in  the  limit  u(h)L  «  1,  for  which 
Ae(x)  ~  B(h)  Ah[— - — ^-L - 1  . 

This  agrees  with  (48)  only  if  W  >  d  -  the  line  width  is  greater  than  the  spacing  - 
or  [u(h)l/wJ  «  1  -  the  peak  of  the  absorption  coefficient  at  the  line  center  is 
small  compared  to  L-1.  If,  on  the  other  hand,  W  «  d  ,  we  obtain 

ft 

AeOO  —  B(k)  Ah  /tt(|-)  JT(k)L  , 

00  v  h 

corresponding  to  the  isolated  line  Laderiburg  -  Reiche  approximation  (Ref.  1,  Eq.  39)* 
For  the  Og  Schumann-Runge  system  at  normal  density  and  temperatures  6000- 
10,000°K  we  have  computed  the  self -absorption  corrections  (Ac(k)/Ac  (h))  of 
Eq.  47  and  48  above.  Purely  Lorentz  broadening  was  assumed,  and  we  have  used 
the  line  widths  of  Breene  (Ref.  5^,  1968,  Vol.  II).  For  path  lengths  L  =  1.0  cm, 
the  maximum  reduction  in  the  band  emissivity  is  found  to  be  ten  percent  at  6000°K 
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and  less  than  one  percent  at  8000°K.  The  results  presented  here  ignore  self¬ 
absorption  effects  in  the  molecular  bands. 


XVE.  Molecular  Absorption  Computer  Programs 

For  the  purpose  of  generating  and  analyzing  the  molecular  absorption 
coefficients  presented  herein,  three  separate  computer  programs  have  been 
developed. 


1.  The  computer  program  MORSE  evaluates  Franck-Condon  factors  q^.,  v„  and"r- 
centroids"  rvty„  (Eq.  36),  assuming  MORSE  potentials.  These  results  have  been 
used  to  extend  existing  tabulations  and  the  extensive  tables  of  RKR-values  kindly 
provided  us  by  Capt.  John  Generosa  of  the  Air  Force  /eapons  Laboratory  (Ref.  37). 

The  formulation  of  Nicholl3  is  closely  followed  (Ref.  43).  For  large 
values  of  the  vibrational  quantum  number,  however,  increased  numerical  accuracy 
is  achieved  by  replacing  his  recursion  relation  for  Laguerre  polynomials  (Ref.  37, 
Eq.  ll)  by  a  relation  we  derived  for  the  ratio  cf  polynomials  of  successive  order. 
With  Lv(Z)  =  (Z),  the  associated  Laguerre  polynomials,  we  obtain 


v-1 

=  2  -  (K~2v)  £  T„,  i  =  1. 
r=0  *  u 


where  Tft  =  1, 


T 


v-r 


Employing  these  relations,  the  code  MORSE  obtains  results  agreeing  with  existing 
MORSE  potential  tabulations  (Ref.  41,  36,  43).  (The  Franck-Condon  factors  snd 
r- centroids  are  punched  on  cards  suitable  for  input  to  code  ABSMOL.) 

2.  The  computer  program  ABSMOL  evaluates  the  averaged  absorption  coefficient 
of  Eq.  40  for  a  particular  band  system  in  both  the  "PQR-Branch"  and  "Smeared-Band" 
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models.  The  temperature  and  spectral  constants  of  Tables  III  and  IV  are  entered 
from  punched  cards,  as  are  the  Franck-Condon  factors  and  r-centroids.  If 
employed,  we  also  enter  parameters  defining  the  dependence  of  the  electronic 
. ransition  moment  Rg  (rv»v, ■)  on  the  r-centroid  ryl v„ . 

The  band-head  energy  and  f-number  for  each  vibrational  transition  and  the 
spectral  absorption  coefficient  and  cummulative  Planck-Mean  generated  by  both 
band  models  are  evaluated  and  listed.  The  histogram  of  line  strength  for  the 
"?QR-Branch"  model  is  accumulated  and  output  on  punched  cards,  and  the  corres¬ 
ponding  absorption  coefficient  is  output  both  graphically  and  on  punched  cards 
suitable  for  input  into  RATRAC  or  COMBDS. 

3.  The  COMbined  BanDS  computer  program,  COMBDS,  was  developed  to  compress 
the  single  band  output  of  ABSMOL  into  an  absorption  coefficient  and  histogram 
for  a  single  molecular  species.  The  output  are  both  punched  cards  suitable 
for  RATRAC  and  graphic j  the  molecular  absorption  coefficients  of  Fig.  4-38 
were  obtained  from  this  code. 

Although  the  combining  of  the  contributing  absorption  coefficient  is 
trivial,  combining  the  separate  histograms  is  complicated.  This  is  because 
each  histogram  describes  the  line  strength  distribution  a  single  band 
in  each  spectral  interval  relative  to  th^  mean  st*  -  ;th  in  hat  interval 
for  that  band.  Combining  sever il  bands  of  differing  mean  strength  into  a 
distribution  about  the  resulting  mean  strength  is  a  non-trivial  exercise  in 
computing  and  it  is  solved  in  the  present  program. 
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V. 


Results  and  Conclusions 


A.  Molecular  Absorption  Coefficient 


The  results  of  our  studies  are  presented  in  both  graphical  and  tabular  form. 


The  molecular  bands  we  have  treated  and  their  spectroscopic  constants  are  in  Tables 
III  and  IV.  Table  V  contains  the  Planck  Mean  per  molecule  (in  units  10  cm  )  for 
each  individual  band  system  (Eq.  5)  -  a  function  of  temperature  only.  Mo.'ecular 


continue  are  not  included  in  this  table. 


To  indicate  the  relative  contributions  of  individual  species  to  air  radiation 


we  present  in  Tables  VI  -  X  the  total  Planck  Mean  -  lines  and  continue  for  temper¬ 
atures  6000  -  24,000°K,  densities  normal  to  10~^  normal  and  broken  down  by  species. 

The  molecular  absorption  coefficients  generated  according  to  our  "PQR-Branch" 
model  of  Chapter  IV  are  in  Figures  4-38.  The  quantity  graphed  is  an  effective 
cross  section  (absorption  coefficient  per  molecule)  in  units  (10~  cm  ),  averaged 
according  to  Eq.  40  over  the  energy  interval  tn  of  Table  III.  These  cross  sections 


are  functions  of  temperature  only.  These  graphs  were  generated  by  the  computer 
program  COMBDS  (Section  IVE)  and  combine  the  individual  band  results  obtained 
by  ABSMOL,  (Similar  graphs  for  the  individual  band  contributions  are  available 
but  are  too  numerous  for  inclusion  in  this  report. ) 


B.  Negative  Atomic  Ions 

To  the  molecular  absorption  coefficients  illustrated  we  have  added  the 
absorption  due  to  the  atomic  negative  ions  0”  and  N”.  The  0~  cross  section  is 
from  Churchill,  et  al.  (Ref.  51) ,  and  the  N~  cross  section  is  from  Kobzev  (Ref.  57) 
an  absorption  process  of  considerable  controversy.  Avilova,  et  al.  (Ref.  4)  use 
an  N~  cross  section  an  order  of  magnitude  greater  than  we  employ  here.  (We 
followed  Avilova's  results  in  our  earlier  report  (Ref.  7).)  Some  recent 
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Russian  experiment'*  1  work  (Ref.  57)  obtain  values  approximately  twice  those 
we  employ  here.  (The  latter  reference  reviews  most,  earlier  work.) 

In  all  cases  the  experimental  procedure  is  one  of  subtraction.  The 
emission  intensity  is  observed  in  spectral  intervals  between  listed  spectral 
lines,  compared  with  a  synthesised  intensity,  and  the  difference  is  attributed 
to  N”.  The  variation  of  the  excess  radiation  with  thermodynamic  variables  is 
consistent  with  a  binding  energy  of  ~  1.0  eV;  only  the  large  magnitude  of  the 
resulting  cross  section  has  been  difficult  to  explain. 

As  we  have  several  times  emphasized  in  this  report,  a  synthesis  of  the 
absorption  coefficient  where  it  is  large  seems  quite  reliable,  but  very  question¬ 
able  where  it  is  small.  As  noted  in  Ref.  57  the  experimentally  accessible 
region  is  rich  in  spectral  lines.  The  line  atlases  we  employ  contain  many  more 
lines  in  this  region  than  do  most  line  locators  -  mostly  lines  of  large  principal 
quantum  number  and  hence  considerable  width.  Asinovskii,  et  al.  (Ref.  57  )  demon¬ 
strate  the  sensitivity  of  their  results  to  the  shift  of  nearby  photoionization 
edges  -  an  effect  roughly  accounting  for  the  presence  of  these  broad  lines. 

Furthermore,  all  the  so-employed  synthetic  spectra  with  which  we  are 
familiar  approximate  the  line  emission  by  an  isolated  line  approximation  (Ref.  l). 
Such  treatments  can  in  no  way  approximate  emission  from  the  wings  of  lines,  which 
we  find  may  be  of  considerable  importance  at  the  temperatures  and  densities  of 
interest.  Finally,  free-free  on  neutral  species  and  molecular  contributions 
mus-c  be  more  carefully  assessed.  For  these  reasons,  we  are  not  yet  ready  to 
accept  the  magnitude  of  the  experimental  N"  cross  section  (Ref.  57  )• 

The  terra  diagram  of  N~  is  in  Figure  53.  Schaeffer  and  Harris  (Ref.  58  ) 
have  estimated  the  ”4)  state  to  be  bound  by  0.844  eV  and  the  *\s  state  to  be 
unbounded.  We  shift  Kobzev* s  estimate  (Ref.  57  ,  Fig.  12,  Curve  (2))  to  this 
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4 

binding  energy  and  discard  his  S  cross  section.  When  we  have  included  the  line 
wing  formulation  of  this  report  in  our  computer  program  RATRAC  we  hope  to  study 
these  questions  further. 

C.  Emlssivity  Results 

The  emissivity  results  obtained  from  our  RATRAC  code  are  presented  in 

n  -2 

Figures  39-46  for  path  lengths  10  -  1.0  cm,  temperatures  0.5  to  5.0  eV  and 

-4 

densities  1.0  to  10  normal.  The  molecular  contribution  Js  included  as  described 
in  Section  IVD,  above.  Finally,  Figures  49-51  present  the  emissivity  at  a  path 
length  L  =  1.0  cm  and  pressures  1.0,  0.5>  and  0.1  atmospheres  as  a  function  of 
temperature. 
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APPENDIX  A 


APPLICATION  OF  SECOND  QUANTIZATION  TO  RACAK  ALGEBRA 
OF  EQUIVALENT  PARTICLES 


AI.  Introduction 

The  use  of  second  quantization  in  atomic  spectroscopy  has  been  in¬ 
creasing  rapidly  in  recent  years.  In  deriving  new  relations,  and  rederiving 
old  ones,  this  method  is  often  more  straightforward  and  more  powerful  in 
taking  account  of  the  PauLl  exclusion  principle  than  working  with  an  anti- 
symmetrized  configuration  space  wave  function  approach.  Consequently,  it 
is  hoped  that  the  explicit  evaluation  of  off-diagonal  matrix  elements  of 
multi-electron  atoms  (necessary,  for  example,  for  calculating  the  auto- 
ionization  cross  section)  is  more  quickly  handled  by  a  computer  if  the 
programing  is  based  on  second  quantization  rather  than  on  fractional 
parentage  coefficients  (FPC)  and  configuration  space  formulae. 

An  exposition  of  the  most  pertinent  formulae  is  given  below,  and  in 
the  last  section  these  are  applied  to  a  matrix  element  whose  value  has  been 
in  dispute. 

All.  Second  Quantization  Formalism  for  Fermions 

This  section  will  begin  with  a  brief  review  of  equivalence  of  the  second 
quantization  (SQ)  method  and  the  configuration  space  (CS)  approach  in  atomic 
problems  where  the  number  of  particles  is  a  constant  of  the  motion.  3he  origins 


l»  n 

— T  {  ~ 


of  SQ  method  will  not  be  described  here,  but  they  can  be  found  in  a  sizeable 
number  of  textbooks  on  the  theory  of  quantized  wave  fields  (Ref.  l). 

The  CS  method  uses,  in  practice,  sets  of  single-particle  states,  and 
the  SQ  method  requires  them  even  in  principle.  In  this  connection,  it  is 
crucial  that  the  Rauli  Principle  is  obeyed. 

The  CS  state  function  of  the  "particles"  generally  consists  of  a 


superposition  of  determinantal  wave  functions 


<Pl(*l)  91(^3)  •••  ‘PqCV 

Y(X]x gX3...xN)  =  x  cp2(x1)  q?2(x2)  ... 


(A.1) 


where  the  ^(x),  cp2(x)  ...  are  a  complete  orthonormal  set  of  Pauli  eigen 
functions  of  a  suitable  single-particle  problem.  For  applications  to  atomic 
spectroscopy  the  fcp}  are  the  solutions  to  the  wave  equation  with  a  screened 
central  field  coulomb  potential  (possibly  non-local  due  to  exchange),  i.e. 


*0^  *  8n*(r>  W*’ 

s 

where  is  the  radial  wave  function,  Y^m  is  a  spherical  harmonic,  and  xm  is 
a  Pauli  spinor. 

In  the  SQ  method  the  state  function  carries  only  occupation  number 
information  regarding  the  states  ^(x),  <f2(x)  ..„  .  For  particles  obeying 
the  Pauli  Principle  (Fermions)  the  occupation  number  of  any  non-degenerate, 
single-particle  state  a  cannot  exceet  =  1.  Thus,  there  exists  a  creation 
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operator  such  that 


B.+  \ct>=° 

i  °  >  =  I  <*  > 


(A. 2) 


where  |  0  >  and  |  a  >  denote,  respectively,  the  states  with  NQ  =  0  and  1.  In 

order  to  assure  the  correspondence  of  the  state  j  Qt,p,  ...  Xu  ...>  = 

a+  at  ...a+  a+  ...  |  0  >  of  SQ  with  the  CS  determinental  state  (A.l),  which  is 

a  0  Yu 

antisynmetric  under  exchange  of  states,  it  is  furthermore  required  that 


a*  +  a*  -  0  (A. 3) 

The  Hermitian  conjugate  of  a+  is  the  annihilation  operator*  with  the 
properties 

aa  I  P  >  =  0  ifpyEa  (A.4) 

aa|o>=|0> 

and 

aa  ap  +  ap  aa  =  0  (A*5) 

The  equations  (A. 2)  and  (A.4)  therefore  imply 

aa  aa  I  0  >  '  0 

4aala>  =  la> 

which  makes  a*  aQ  the  occupation  number  operator  for  state  a. 


.I»Q. 


Moreover,  in  order  to  maintain  the  equality  of 


i  1  f  5  5 

'Pa<x2> 

<f7(*l)  <PY(xg) 

<ap|Y6  >  =  ?J  ‘S*3* 

?B(xl)  <Pg(x2) 

<P5<X1) 

6ory  6ps  ’  6a6  6pY 

and 

<  a  p  |  y  6  >  -  <  0  |  ap  aQ  a+  a+  |  0  >, 

it  is  required  that 

aa  ap  +  ap  aa  =  6ap  (A.  6) 


Next,  operators  and  their  matrix  elements  will  be  compared  in  the  CS 
and  SQ  approaches.  If  the  "particles"  are  identical,  operators  must  be  sym¬ 
metric  under  exchange  of  particle  labels  in  the  CS  representation.  Operators 
considered  are  made  up  of  a  sum  of  individual  pa-tiele  operators  of  the  same 
form, 

F  =  H  f<xi) 

i 

and  operators  consisting  of  a  sum  of  pair  contributions  of  the  same  form 


0  =£  e(xr 

ij 


If  operators  of  the  F  and  G  type  are  used  to  compute  CS  matrix  elements 
with  fully-antisymnetrized  wave  functions  (A. l), 
found  in  SQ  by  using 


f  A 


then  the  same  values  will  be 


p  |  f  1  a  >  aff 


2$=£  apa^<P<JS8,.TU>auaT 

U  POTU 


(A.  7) 


'f'  *4*  | 

and  the  corresponding  normalized  state  vector  ag  ...  a^  j  0  >  . 

In  (A. 7)  the  matrix  elements  are  defined  6s 

<  P  1  f  |  cr>s  J d3x  <p*(x)  f(x)  cpCT(x) 
and 

Jo  O  *  ft 

d^Xg  p^Xj)  q>ff(xg)  gi^Xp)  ^(xj)  q^Xg) 

The  sums  p,o/r  ,u  extend  over  a  complete  set  of  single-particle  states.  Thus,  in 
SQ  the  Pauli  Principle  is  satisfied  hy  using  the  anti-commutation  properties  of 
the  a+  and  a  operators. 

It  is  worth  noting  that  the  equivalence  of  the  two  approaches  is  formally 
rooted  in  the  correspondence  of  the  Y  of  (A.l)  and 

t(x) 

ft 

In  other  words  if  the  7  and  Y  in  CS  matrix  elements  are  replaced  hy  the  operators 
f(x)  ans  if  (x),  and  the  matrix  elements  vith  respect  to  the  occupation  number 
state  vectors  are  taken,  the  result  will  be  the  same  as  the  CS  calculation. 
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AIII.  Rotational  Transformation  Properties 

In  order  to  be  able  to  exploit  rotational  symmetry,  a  discussion  ol’  it 
in  the  context  of  second  quantization  will  be  given  in  this  section. 

Consider  two  coordinate  systems,  "prime"  and  "unprime"  related  to  each 
other  through  a  rotation  specified  by  a  set  of  parameters  R  such  as  the  three 
Euler  angles.  In  order  that  the  form  of  the  equations  of  motion  and  the  arrti- 
comnutation  relations  he  preserved  in  the  two  coordinate  system,  there  must 
exist  a  unitary  operator  U  such  that 

*  -f  I-  +  t»“l 

a  =  U  a  U 

K  K 

< 

Another  way  of  relating  a  and  a  is  to  examine  the  commutation  relations 

K  H 

of  a+  and  the  orbital  angular  momentum  operator 

Xj  =  J^a*  <  p  |  Lj  j  a  >  aff>  (A.8) 

Let  the  greek  subscripts  denote  the  principal  quantum  number,  n,  total  orbital 
angular  momentum,  l,  total  spin  angular  mane  return  l/2,  z  component  of  the  orbital 
angular  momentum,  m^,  and  z  component  of  the  spin  angular  momentum,  ms.  Then, 
using  (A. 3),  (A. 5),  (A. 6),  and  (A.8),  we  obtain  commutation  relations 


r  »  + 

CV  =  H  aY 

C X+,  =  [X(X+1)  -  nfc(mxil)3l/2a^  , 


(A.9) 


where  £+  =  (X^  £  iXy)//2  and  Y+  =  {n,  X,  ra^  +  1,  ms}.  Thus,  the  commutation 

relations  of  X.  and  a+  are  the  same  as  those  of  L.  and  Y.  .  Since  the 
J  3  Jim 
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«  + 

commutation  relations  generate  the  local  values  of  a  ve  know  their  values  in 
the  large.  Just  as  with  the  Ylm‘  For  the  latter 

'  £  bL.(>»  W«) 

vhere  D^mt (R)  is  the  matrix  element  of  the  1-dimensional  irreducible  represent¬ 
ation  of  the  rotation  R  for  the  three-dimensional  rotation  group.  If  spin  is 
now  considered  by  Introducing  the  spin  angular  momentum  operator 


fc?J  ap  <  Pi  SjIt  >aT 

and  its  commutation  relations  with  a+,  one  readily  sees  that  a+  transforms 
according  to  a  s  =  l/2  representation  in  its  spin  quantum  numbers.  Therefore, 


DvH‘1=ZDWB)I>"/»’<R)a“ 

?  i  8  S 


m  m 
e  s 


Because 


-  V8'1’ 


the  relation  complex  conjugate  to  (A.10)  becomes 


(A.  10) 


o  V”'1  Bm'm  f8'1)  <A-U> 

11  S  s 

In  other  words,  aQ  does  not  transform  according  to  l/,a(R)  E^/^(R). 

Relations  (A.10)  and  (A.ll)  are  useful  in  establishing  the  transformation 
properties  of  more  complicated  quantities  built  up  from  sets  of  products  of  a+ 
and/or  a,  as,  for  example,  multi-particle  creation  operators. 
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If  the  spin-orbit  interaction  is  neglected,  as  will  be  the  case  through¬ 
out  this  appendix,  each  stationary  state  of  a  system  of  bound  atomic  electrons 
will  have  an  integral  value  of  the  total  orbital  angular  momentum,  L,  and  a 
half- integral  or  integral  value  of  the  total  spin  angular  momentum  S.  We 
show  how  to  combine  products  of  the  a+l s  to  form  such  states. 

In  combining  two  angular  momenta  and  Jg  so  as  to  form  an  eigenstate 
of  their  total  angular  momentum  J  (with  z  components,  m^,  mg,  and  M,  respect¬ 
ively),  the  vector  coupling  coefficient  <  m^  Jg  |  J  M  >  appears  according 
to 

1  ^  >  a^<  ift  3ft  1  JM  >  I  Jl®!  >  I  J2®2  >  (A.  12) 

“l^ 

The  explicit  form  of  the  vector  coupling  coefficient  is  given  in  numerous  books, 
e.g. ,  Brink  and  Satchler  (Ref.  l),  and  Judd  (Ref.  4).  However,  for  present 
purposes  the  following  relations  (Ref.  l)  will  suffice. 

y]  <  «)g®g  |  JM  >  <  JM  ]  Jgn^  >  =  6(11^,11^)  6(mg,n^) 

JM 

(A.  13) 

JM  |  3ft  j2mg  >  <  3ft  Jgmg  |  J*M«  >  =  5(^P  )5(J,J* ) 

c  3ft  3ft  !  JM  >  =  (-1)J1+VJ  <  jgmg  3ft  I  JM  >  (A.  14) 

If  and  jg  refer  to  the  angular  momentum  of  identical  Fermions,  then  (A. 12) 
must  be  brought  into  accord  with  the  Pauli  Principle.  This  is  done  by 
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constructing  states 


I  JXJ2  JM  >a  -53<3l"l  Vfe  1  M  >  a3j^  1  °  > 

“l"? 


(A.  15) 


s0Si°> 

xo  prove  that  the  state  |  J^Jg  JM  >a  actually  has  the  eigen  value  J  we 
examine  the  behavior  of  0^  under  rotation.  The  rotational  behavior  of  tne 
creation  operators  is  given  by  (A.10),  written  explicitly  as 


S.'""1  >E  Di'<s)  % 


Hence, 


°  <4 u'1  =Eff  V)/1  <  3  a  3&‘\;k> 

=?  4  4 <  j * 3*  i  ■ » >  x 

"iW 


(A.l6) 


(A.  17) 


J,  d 


1>n1m1  Dngmg  Jlml  J2“2  ^  JW'  >  DNfN  <  **  l  Jlml  >  (A.lB) 


(Ref.  1). 


Using  the  orthogonality  relation  (A. .13}  and  (A.18),  Eq.  (A. 17)  becomes 

U  <4  U‘1  L"  31°1  Vte  1  JM'  >  X, 

E  <6  dm'« 
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:  **  Jv»+i4  '■'Vi*  •  ■ 


Thus ,  c£  transforms  accorOinj  to  the  J-dimensiona'..  irreducible  representation 
of  the  rotation  group. 

From  eq.  (A. 15)  on  we  have  treated  only  angular  momenta  and  Jg  of 
two  identical  particles.  In  the  present  approach,  where  L-S  coupling  is 
ignored,  each  single-particle  state  is  an  eigenstate  of  spin  and  orbital 
angular  momentum  (total  and  z  components).  This  causes  only  a  minor  compli¬ 
cation.  Since  spin  and  orbital  motion  are  uncoupled,  they  transform  separately, 
according  to  (A. 19).  The  relation  equivalent  to  (A.19)  becomes 


m’t' 


where 


-LS 
°mt 


<  im£,  jt’mJ  |  IM  >  <  1/2  m£ 


1  m. 

2  ms 


|  St  > 


+  + 
a  a  . 
v  v' 


(A.20) 


States  of  more  than  two  particles  having  a  certain  total  L  and  S  can  thus 
be  built  up  from  pairs  of  particles  with  given  transformation  properties  (and 
a  single  particle)  if  the  number  of  particles  is  even  (odd).  For  example,  a 
four-electron  state  of  particular  L  and  S  can  be  constructed  from  pairs  of  two- 
electron  eigenstates  as  follows 

1  Vl  *2S2  (L’s‘)>  X3S3  V4  <L"S">i  ItS  >  - 

J]  <  L»M*  L"M"  I  LM  >  <  S'T'  S"T"  |  ST  >  0^'  0^"  |  0  > 

with  the  operators  given  by  (A.20);  some  of  the  z  component  labels  ,  as  well 
as  the  normalization  constant,  have  been  suppressed. 


For  the  presentation  of  elementary  angular  momentum  coupling  relations  > 
such  as  those  above,  the  vector  coupling  coefficient  notation  suffices.  However, 
for  problems  involving  a  great  many  angular  momenta  (as  in  the  following  section), 
it  is  simpler  to  exploit  the  symmetries  of  the  3-J  symbol,  which  is  defined  by 
the  relation 

<  a  a  b  p  !  c-v  >  (-l)a“b"Y  =  /2c+I b  A  (A. 21) 

As  a  result  of  (A. 21)  the  3-J  symbol  has  the  following  very  useful  pro¬ 
perties.  For  non-cyclical  permutations  of  columns, 


=  (-1)J1+J2+J3 


/h^h  \=(.i)31+32+35  /Ji3233\ 
VWV  VWj  / 


(A. 22) 


(A.23) 


Numerous  authors  have  introduced  their  own  transformation  symbol  analogous  to 
(A. 21).  As  a  result,  errors  in  phase,  due  to  different  conventions,  can  easily 
be  made  by  the  beginner.  A  good  example  of  this  is  the  different  uses  of  the 
same  symbol  V  by  Racah  (Ref.  5),  and  de  Shalit  and  Talmi  (Ref.  7).  Probably 
the  most  complete  listing  of  the  various  transformation  symbols  is  collected 
in  Yutsis,  et  al.  (Ref.  8).  The  3-J  symbol,  however,  seems  to  be  becoming  the 
most  widely  used. 
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which  in  turn  is  substituted  into  (A, 25),  and  one  obtains 


<  a  b  j  V  |  c  d  >  <Va>  Vc>°k  <4aV  Vb>  * 

31  ^Va'  Vb'  Vc>  "aV1(\“b'  VV  * 


(A.  26) 


x  6(ms(a),  bb(y))  6(afl(p), 


Definitions  introduced  in  (A. 25)  are  the  following: 

CK(X  m  l'm,)  =  sind  8^0^,  9Vmt 


where 


and 


Y  *  -i-elnx P0. 
Jtm  /2,x  im 


”  •  JC 


^Va'Vb-Vc'Va)  “  eY /  t4rT\W  V8>  V1)  *4<*)^x** 

•O^o  rx+1 


with  the  R  subscript  e  =  {n  i  } . 

It  can  be  shown  that 

=  ( -l)m' ( ( 2A+l)(2t'  +1)] (J  4'  fj  J.B)  (A. 27) 

and  thus  the  angular  integrations  are  disposed  of.  The  evaluation  of  will  not 
be  discussed  further  because  the  evaluation  of  these  integrals  is  not  related  to 
the  influence  of  the  Pauli  Principle  on  the  rotational  properties  of  matrix  elements. 
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The  task  of  finding  the  non-vanishing  parts  of  the  various  products  of 
annihilation  and  creation  operators  in  (A. 24)  can  he  systematically  carried 
out  with  Wick’ s  Theorem  (Ref.  6),  or  by  matching  directly  creation  and  an¬ 
nihilation  operator  pairs  so  that  there  are  no  non-vanishing  matrix  elements. 
It  seems  that  this  procedure  is  not  hard  to  do  on  a  computer. 

To  some  extent,  this  investigation  arose  from  .1  discrepancy  in  the 
value  of 

X  =  <  2s2  (1S)  2p2  (3P)j  3P  |  V  |  2s3d(1D),  2p2(3P);  3P  > 

calculated  in  two  attempts  to  apply  the  CS  approach  based  on  FPC,  A  computer 
program  gave  X  =  X+,  where 

X+  =  i  R1  (2s,  2p,  2p,  3d)  +  §  R2(2s,  2p,  3d,  2p) 


while  an  application  of  Fsno*  s  prescription  (Ref.  3)  gave  X=X-.  It  was  decided  to 
evaluate  X  by  an  independent  method,  i.e.  by  the  use  of  SQ.  The  main  steps 


of  this  approach  will  be  sketched  next. 

The  eigenvalue  of  the  z  component  of  the  orbital  and  spin  angular 
momenta  we  denoted  by  cr  =  0  and  for  an  s  electron,  tt  and  for  a  p  electron, 
and  A  and  for  a  d  electron. 

Due  to  (A. 15)  and  (A. 21),  the  normalized  ket  is  given  by 


lov. 

~  fo  2-i  cmrr’A  crrm' 


+  +  ! 
a^  a.  a. 


where,  for  simplicity,  the  spin  <7  labels  have  been  omitted  on  the  creation  oper- 

/  g  }  f  g  } 

ators  as  well  as  on  A'  '  and  B'  7 . 


-60- 


_(L)  ,  ,  /0  2  2\  /ll  1\/2  1  1  \ 

OTTTT'A  01  5  yO  A  -A  f  \TT  TT*  -(n +n*  )  J  \A  TT+tt*  -(AtHWC'  )y 

/  .A-3  A  s  °  Vi  i  1  V°  1  1 

Y°a  aA  “(aa+  aA^ y  yarr  arr'  “(aa  +  atT»  y  c^+a^,  "(an+  cr^ 


Similarly,  the  normalized  bra  is  given  by 

<  A  I  =|  £<  0  I  aaa  W*  Aan''TTH‘  A^W 


with 


,S 

on"TT=* 


/i  i  °VI  i  1  V°  1  1  ) 

\°  °/\  V  CTtt  “(ffTt"+  aTT*  ^  '  0  V,+  aTT*  ”K"+  0TT'"  ^ / 


/l  1  1  \  /0  1  1  \ 

AanV»  ^tt"  u«  -(Ti"-Kr- )/  \0  7T"+n- 


and 


a*  <  a  p  |  V  I  Y  5  >  a§  ay 


Non-vanishing  contributions  of  V  to  X  occur  if  the  two  creation  operators 
involve  an  s  and  p  electron,  and  if  the  two  annihilation  operators  refer  to  d  and 
p  electrons.  The  contribution  to  the  electron-electron  coulomb  energy  matrix 
element  X  vanishes  if  particles  other  than  p  electrons  are  both  "created  and 
annihilated"  in  V.  Thus,  it  turns  out  that  the  total  number  of  non-vanishing 
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<  A  |  a*  a^  afi  j  B  >  elements  la  sixteen. 

After  carrying  out  all  contractions  of  a+  and  a,  X  may  be  written  in  the 

form 

X  =  23/2  (II  -  I) 

where 


n‘ZAX-  AX*  <”•  •civi"’  ^ 

1  *2  A0*n'  Aoftn>  <  ■°  I  v  i  A>  "  >  B0inni 


(A.  31) 


Only  a  single  term  in  both  <  n,  -o  |  V  |  n,  A  >  and  <  tF,  -a  |  V  |  A,  n  > 
expansions  does  not  vanish  (see  A. 26),  as  can  be  seen  from  the  explicit  values  of 
the  3-j  symbols  (Ref.  l)  occurring  in  (A  )  and  (A. 27),  thus 


<  ti,  -o  |  V  |  A,  tt  >  =  |  (-l)n  6(0^.,  a^)  6(0^^)  6(ff,TT4A)Ri 

<  ",-c  |  V  !  tt,  A  >  (-1)^  *5(o^;  an)  6 (o-q>  aA)  6(ff,n+A)R2 


A  further  simplification  takes  place  at  this  stage  because  X  is  diagonal 
In  the  z  components  of  the  total  spin  and  orbit  angular  momentum  -  in  fact,  X 
is  independent  of  both.  Therefore,  both  of  these  conqponents  may  be  taken  equal 
to  zero. 

Then  one  finally  finds  that  the  sums  of  (A.3l)  give 


I  =  -rV6  /2 

Il/l  =  6/5 ‘R2/!*1 


■62- 


Therefore, 

X  *  X_  =  g  R^s,  2p,  2p  3d)  -  |  R2(2s,  2p,  3d,  2p) 
■which  is  the  result  obtained  by  Fano*  s  treatment.* 


* - - 

Dr.  Paul  Bagus,  Private  Communication,  IBM  Research  Laboratory,  San  Jose,  Calif. 
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AP3FENDIX  B 


ON  THE  NUMERICAL  INTEGRATION  OVER  DETAILED  SPECTRAL  LINE  IROFILES 

The  preliminary  emissivity  results  reported  in  Ref.  1  have  been  found  to  b 
numerically  unstable  for  varying  choices  of  the  photon  energy  mesh  interval 
at  densities  below  10  normal.  At  such  low  densities  collisional  line 
widths  may  become  less  than  10  eV  and  any  computationally  practical  mesh 
interval  Aw  may  contain  a  large,  but  unknown,  number  of  strong  narrow  lines. 

To  treat  these  in  the  evaluation  of  the  contribution  to  Eq.  4  from  the 
photon  interval  [w,  w  +  Am],  we  write 

p’(u>)  *=  f(w)  +  g(u>) 

where  f(u>)  consists  of  all  of  the  contributions  to  p’fco)  which  vary  slowly 
in  the  neighborhood  of  <u. 

Then 

Ac  (w)  =  -  e"4 '  (^ ' )  h-j 

0)  oT 

«  tm  [1  -  e”f(u,>L  -A  f  %uVe(tt},)L} 

oT  w 

As  described  in  Section  IV  E  of  Ref.  1,  f(w)  includes  continuous 
absorption,  "broad"  lines  (width  W>  Aw),  "weak"  lines  (absorption  at  line 
center  less  than  continuum  absorption) ,  and  the  wings  of  "strong  narrow" 
lines  (with  line  centers  lying  outside  the  interval  [»  -  Aw,  w  +  2  Aw]) . 
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We  visualize  the  situation  as  in  Fig.  1.  The  lines  numbered  (l)  through 
(6)  are  all  "strong  narrow"  lines  by  the  a’jove  criteria.  The  integration 
over  the  interval  [to,  w  +  Aui]  is  performed  by  subdividing  Aw  into  unequal 
subintervals  sufficient  to  resolve  in  [w,  w  +  Aw]  all  the  "strong  narrow" 
lines  whose  centers  lie  in  the  larger  interval -  Aw,  w  +  2  Aw],  accumulating  at  this 

submesh  the  total  absorption  coefficient  due  to  all  lines  in  the  larger 
interval  (plus  f(w))  and  integrating  by  trapezoidal  integration. 

The  earlier  numerical  instability  referred  to  above  was  due  to  an 
unsuitable  prescription  for  defining  the  energy  submesh.  With  extensive 
computer  experimentation  we  have  obtained  a  prescription  which  is  efficient 
and  accurate  for  all  cases  we  have  tested  (densities  down  to  10”^  normal 
and  temperatures  0.5  eV  to  100  eV)  . 

Two  passes  arc  made  through  the  contributing  strong  narrow  lines. 

Every  such  line  whose  absorption  coefficient  at  line  center  lies  below  the 
wing  of  any  other  line  is  removed  from  the  mesh-defining  considerations. 

In  the  example  of  rig.  1,  only  lines  (1) ,  (3) ,  and  (4)  will  define  the 
submesh.  The  minimum  in  the  absorption  coefficient  due  solely  to  adjacent 
pairs  of  these  remaining  dominant  lines  is  evaluated,  assuming  Lorentz  wing 
profiles  -  resulting  in  the  energy  values  xg  and  in  the  figure.  The 
width  of  each  dominant  line  then  is  employed  to  subdivide  the  mesh  between 
its  two  adjacent  minima  -  line  (1)  governs  the  interval  [xQ  =  o>,  xg]  and 
line  (3)  the  interval  [xg,  x1  =  w  +  Aw].  Specifically,  the  interval 
[x2,x.]  is  subdivided  by  the  energies  w-  ±  B  .  where  B  =  0.25  W,  and 
Bn  =  1:6  Bn^1(n  =  1,2, . . .) ,  u>3  being  the  line  center  of  line  (3)  and  W^ 
its  width.  The  energy  submesh  so  obtained  is  then  ordered,  the  absorption 
coefficient  due  to  all  lines  evaluated  over  it,  and  the  integration  is 
performed. 
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The  submesh  defined  in  this  fashion  by  definition  resolves  the 
dominant  lines  and  by  experience  is  sufficiently  fine  grained  to  resolve 
the  non-dominant  ones.  Results  included  in  this  report  were  evaluated  by 
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Table  II.  SOME  Nil  LINE  WIDTHS 


17  with  Bely- Griem- Coulomb  Gaunt  factor  (Fig.  5) 
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Table  IV.  MOLECULAR  SPECTROSCOPIC  CONSTANTS 
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Table  VI.  TOTAL  PLANCK  MEAN  (CM-1)*-  (By  Species) 
Density  p/pQ  *  1.0 


TEMPERATURE  (ft) 


SPECIES 

4000 

6000 

8000 

10.000 

12.000 

18.000 

24.000 

*> 

4.14-7 

3.24-4 

1.82-2 

9.92-2 

1.05-1 

8.96-3 

1.69-3 

°2 

6.60-4 

3.61-3 

3.02-3 

2.88-3 

2.69-3 

1.31-4 

1.59-5 

NO 

2.68-4 

9.22-3 

2.80-2 

3.32-2 

1.98-2 

2.09-3 

1.50-4 

< 

2.03-11 

2.26-6 

1.05-3 

1.40-2 

2.38-2 

1.15-2 

2.29-3 

- 

2.43-5 

2.06-6 

- 

- 

• 

- 

N“ 

- 

1.28-6 

1.49-4 

2.32-3 

1.20-2 

- 

•» 

N 

- 

3.26-4 

1.50-1 

3.4l 

1.60+1 

9.55+1 

1.43+2 

N+ 

- 

- 

2.58.6 

1.3-3 

3.73-2 

3.37 

2.78+1 

0“ 

- 

6.17-4 

2.81-3 

1.06-2 

3.06-2 

- 

- 

0 

- 

6.75-4 

2.48-2 

2.04-1 

8.61-1 

7.78 

1.75+1 

0+ 

- 

- 

- 

9.65-4 

1.94-1 

2.65 

TOTAL** 

1.51-2 

2.29.1 

I.78+I 

1.07+2 

1.92+2 

L.  _ 

**  Includes  Free-free  and  higher  degrees  of  ionisation 
*  6.75-4  denotes  6.75xl0~4 


Table  VII.  TOTAL  PLANCK  MEAN  (Clf1)*  -  (By  Species) 

-1 

Density  p/pQ  -  10 


SPECIES 

4000 

6000 

TEMPERATl 

8000 

RE  (fc) 
10.000 

■■■ 

mi 

»s 

4.22-8 

3.11-5 

9.68-4 

1.66-3 

1.15-3 

1.04-4 

1.74-6 

0* 

3.10-5 

4.55-5 

3.55-5 

3.02-5 

2.66-5 

8.21-7 

2.08-8 

NO 

1.79-5 

3.21-4 

6.79-4 

4.58-4 

2.06-4 

1.17-5 

1.77-7 

4 

5.10-12 

1.14-6 

2.63-4 

8.34-4 

8.32-4 

2.15-4 

1.30-5 

no2 

- 

9.50-8 

5.23-9 

-- 

- 

- 

n“ 

- 

7.57-8 

7.36-6 

9.47-5 

3.89-4 

- 

- 

. 

9.98-5 

3.47-2 

4.79-1 

1.76 

6.95 

5.25 

N+ 

- 

2.80-6 

5.79-4 

1.26-2 

9.0-1 

5.16 

0* 

- 

1.32-5 

6.29-5 

3.28-4 

9.48-4 

- 

- 

0 

- 

7.27-5 

2.59-3 

2.82-3 

8.51-2 

6.56-1 

7.32-1 

0+ 

- 

- 

- 

9.14-6 

3.10-4 

5.77-2 

5.68-1 

TOTAL** 

5-94-5 

3-93-2 

4.86-1 

1.86 

8.57 

1.18+1 

**  Includes  Free-free  and  higher  degrees  of  ionization 
*  6*75-4  denotes  6.75x10  • 


Table  VIII.  TOTAL  PLANCK  MEAN  (CM"1)*-  By  Species 

_? 

Density  p/p  Q  »  10 


SPECIES 

4000 

6000 

TEMPERATE 

8000 

RE  (k) 
10.000 

12.000 

18.000 

24.000 

4.31-9 

2.58-6 

2.15-5 

1.89-5 

9.71-6 

6.55-8 

6.67-11 

°2 

6.33-7 

4.84-7 

3.41-7 

2.S7-7 

2.39-7 

1.76-9 

4.65-12 

NO 

8.15-7 

9-51-6 

1.03-5 

4.66-6 

1.81-6 

2.07-8 

3.68-11 

3.50-12 

5.30-7 

2.10-5 

2.82-5 

2.34-5 

1.57-6 

2.14-8 

"°2 

- 

2.91-10 

7.99-12 

- 

- 

- 

- 

N" 

m 

3.87-9 

3.05-7 

5.04-6 

1.08-5 

- 

- 

N 

- 

2.86-5 

5.17-3 

4.71-2 

1.61-1 

2.91-1 

8.35-2 

S* 

- 

- 

1.51-6 

1.89-4 

3.82-3 

1.79-1 

6.35-1 

0“ 

- 

2.42-7 

1.77-6 

1.02-5 

♦ 

2.71-5 

- 

- 

0 

- 

7.47-6 

2.63-4 

2.08-3 

8.01-3 

3.24-2 

1.25-2 

0+ 

- 

- 

1.2-8 

2.89-6 

9.70-5 

1.39-2 

7.60-2 

TOTAL** 

4,95-5 

5.49-3 

4.95-2 

1.73-1 

5.17-1 

8.12-1 

**  Includes  Free-free  and  higher  degrees  ©f  ionisation 
*  6.75-4  denotes  6.75xl0“4 


Table  IX.  TOTAL  PLANCK  MEAN  (CM"1)*-  (By  Species) 

_■* 

Density  p/pQ  »  10 


TEMPERATURE  (ft) 


SPECIES 

4000 

6000 

8000 

10.000 

12.000 

18.000 

24.000 

n2 

4.33-10 

1.42-7 

2.49-7 

1.55-7 

1.73-7 

1.66-11 

7.30-15 

°2 

7.58-9 

4.92-9 

3.38-9 

1.16-9 

4.60-lC 

5.65-13 

5.14-16 

NO 

2.83-8 

2.26-7 

1.10-7 

4.22-8 

1.14-8 

5.90-12 

4,07-15 

N2 

1.89-12 

1.62-7 

7.82-7 

8.17-7 

4.85-17 

3.13-9 

2.30-11 

no2 

- 

6.93-13 

- 

- 

n“ 

•» 

1.63-10 

1.01-8 

8.75-8 

«• 

- 

- 

N 

mm 

6.70-6 

5.55-4 

4.41-3 

1.16-2 

5.11-3 

8.94-4 

N+ 

~ 

5.24-7 

5.80-5 

1.15-3 

2.35-2 

6.63-2 

0" 

- 

4.39-9 

5.46-8 

2.98-7 

- 

- 

- 

0 

- 

7.45-7 

2.60-5 

1.97-4 

7.06-4 

6.42-4 

1.36-4 

0+ 

- 

- 

3.9-9 

9.01-6 

3.10-5 

2.05-3 

7.96-3 

TOTAL** 

7.98-6 

5.83-4 

4.68-3 

1.35-2 

3.13-2 

7.54-2 

**  Includes  Free-free  and  higher  degrees  of  ionization 
*  6.75-4  denotes  6. 75x10” 4 


Table  X.  TOTAL  PLANCK  MEAN  (CM  ■*■)*-  (By  Species) 

_k 

Density  p/p  »  10 


SPECIES 

4000 

6000 

TEMPERATU] 

8000 

N2 

4.37-H 

3.56-9 

2.52-9 

°2 

7.90-11 

4.96-H 

2.04-11 

NO 

9.06-10 

3^9-9 

1.07-9 

4 

1.04-12 

1.69-8 

2.44-8 

no2 

- 

- 

N f 

- 

- 

• 

N 

- 

1.03-6 

5.40-5 

N+ 

- 

3.9-11 

1.63-7 

0” 

- 

- 

- 

0 

7.52-8 

2.71-6 

0+ 

- 

- 

1.31-9 

TOTALS 

1.13-6 

5.69-5 

1.51-9  1  2.89-7  7.66-6  2.21-4  7.5-4 


Includes  Free-free  and  higher  degrees  of  ionization 
6.75-4  denotes  6.75xl0‘4 
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Thermal  Averaged  Electron-Ion  Collision 
Broadening  Gaunt  Factors  G(w,T)(kT  =  2. 


terms!  Averaged  Electron-Ion  Collision  Broadening 
Mint  Factors  (Eq.  18c)  (for  various  temperatures) 
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Fig.  49  Emissivity  vs.  Temperature  at 
Pressure  Balance  (1. 0  atxnoa. } 
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